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Abstract. Black holes are a common feature of the Universe. They are observed as stellar 
mass black holes spread throughout galaxies and as supermassive objects in their centres. Ob- 
servations of stars orbiting close to the centre of our Galaxy provide detailed clear evidence 
for the presence of a 4 million Solar mass black hole. Gas accreting onto distant supermassive 
black holes produces the most luminous persistent sources of radiation observed, outshining 
galaxies as quasars. The energy generated by such displays may even profoundly affect the 
fate of a galaxy. We briefly review the history of black holes and relativistic astrophysics be- 
fore exploring the observational evidence for black holes and reviewing current observations 
including black hole mass and spin. In parallel we outline the general relativistic derivation 
of the physical properties of black holes relevant to observation. Finally we speculate on fu- 
ture observations and touch on black hole thermodynamics and the extraction of energy from 
rotating black holes. 
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1 Introduction 


Black holes are exotic relativistic objects which are common in the Universe. It has 
now been realised that they play a major role in the evolution of galaxies. Accretion 
of matter into them provides the power source for millions of high-energy sources 
spanning the entire electromagnetic spectrum. In this chapter we consider black holes 
from an astrophysical point of view, and highlight their astrophysical roles as well 
as providing details of the General Relativistic phenomena which are vital for their 
understanding. 


To aid the reader in appreciating both aspects, we have provided two tracks through 
the material of this Chapter. Track 1 provides an overview of their astrophysical role 
and of their history within 20th and 21st century astrophysics. Track 2 (in italic text) 
provides the mathematical and physical details of what black holes are, and provides 
derivations of their properties within General Relativity. These two tracks are tied 
together in a way which we hope readers with a variety of astrophysical interests and 
persuasions will find useful. 
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2 A Brief History of Astrophysical Black Holes 


Although the term “black hole” was coined by J.A. Wheeler in 1967, the concept of a 
black hole is over two hundred years old. In 1783, John Michell was considering 
how to measure the mass of a star by the effect of its gravity on the speed of the light 
it emitted. Newton had earlier theorized that light consists of small particles. Michell 
realized that if a star had the same density as the Sun yet was 500 times larger in 
size, then light could not escape from it. The star would thus be invisible. He noted, 
however, that if it was orbited by a luminous star, the measurable motion of that star 
would betray the presence of the invisible one. 

This prescient, but largely forgotten paper, embodies two important concepts. The 
first is that Newtonian light and gravity predicts a minimum radius R = 2GM/¢ for 
a body of mass M from within which the body would not be visible. The second is that 
it can still be detected by its gravitational influence on neighbouring stars. The radius 
is now known as the Schwarzschild radius of General Relativity and is the radius of 
the event horizon of a non-spinning black hole. Black holes are now known to be 
common due to their gravitational effect on nearby stars and gas. Pursuing Newtonian 
black holes further leads to logical inconsistencies and also the problem that relativity 
requires the speed of light to be constant. 

The concept re-emerged after the publication of Eintein’s General Theory of Rel- 
ativity in 1915 when Karl Schwarzschild found a solution for a point mass. Einstein 
himself "had not expected that the exact solution to the problem could be formulated". 
It was not realised at the time that the solution represented an object which would turn 
out to be common in the Universe. Chandrasekhar in 1931 [9] discovered an upper 
limit to the mass of a degenerate star and which implied the formation of black holes 
(although this was not spelled out). Eddington, who wrote the first book of General 
Relativity to appear in English, considered the inevitability of complete gravitational 
collapse to be a reductio ad absurdum of Chandrasekhar’s formula. The concept was 
again ignored for a further two decades, apart from work by Oppenheimer and Sny- 
der who considered the collapse of a homogenous sphere of pressureless gas in 
GR, and found that the sphere becomes cut off from communication with the rest of 
the Universe. In fact, what they had discovered was the inevitability of the formation 
of a black hole when there is no pressure support. 

With this brief introduction to the early history, we now give some details of what 
is now understood by a Schwarzschild black hole. 


2.1 The Schwarzschild Metric 


The simplest type of black hole is described by the Schwarzschild metric. This is a 
vacuum solution of the Einstein field equations in the static, spherically symmetric 
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case, and takes the form 


2GM 1 2GM\ 7! 2 2 sin? 0 
as? = (1 2 ) ae (1 7 ) a? -Gde — "Fae? (21) 
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where M is the mass of the central object. Most of the observational properties of 
black holes that we need follow from this metric, rather than needing the deeper level 
of the field equations themselves to understand them, but to give an indication of the 
achievement of Schwarzschild, and since we will wish to consider briefly more general 
spherically-symmetric black holes as well, we give a brief sketch of how the metric is 
derived. For more details, and a description of the sign conventions we are using, see 
Chapters 8 and 9 of [30]. 
The field equations themselves are 


Gw _ —8r Tw (2.2) 


where the Einstein tensor G is a trace-reversed version of the Ricci tensor Ry, 
and T is the matter stress-energy tensor (SET). If the matter SET is traceless, in 
particular if we are in a vacuum case, then clearly G „ must be traceless as well, and 
so the vacuum equations are that all components of the Ricci tensor are zero: 


Rw =0 (2.3) 


The Ricci tensor is in turn a contraction, Ryy = R> uv Of the Riemann tensor, which 
is defined by 


Ta, Tag 
axl Ox 


where I" „g is the connection. For a general metric gy, defined by 


RY apy = HIH op ya — oI ga (2.4) 


ds? = Juvdz” dx” (2.5) 


the connection is given in terms of the metric by 


CEOE. — | Jay “tas (2.6) 
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We can thus see that solving the apparently simple equation involves poten- 
tially formidably complicated second order partial differential equations, and it is not 
surprising that Einstein considered it unlikely that an exact solution for the metric 
around a spherical body would be found. In the current case, where we have assumed 
everything is static, we can attempt an ansatz for the metric which simplifies things 
considerably. We let 


ds? = Ad? — B dr? — r?° (d6” + sin’ 6 do?) (2.7) 
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where A and B are functions of r alone. This means that our PDEs become ODEs, 
and although the working is still complicated, we eventually find that we require 


A" A’ A’ B' A’ 
Ro = -3p 56 a rao (2:8) 
A" A! A’ B! B' 

Mir = 5A a(S i z) aN (a) 
E E E an 2.10) 
aS B "2B\ A B] ” u 
R3 = Rysin?@=0. (2.11) 


We can eliminate A and find a simple equation for B by forming the combination 


Ro Ru  2Ro 


H H = 2.12 
A B r2 9 ene 
which yields the ODE 
dB B(I-B 
= ( (2.13) 
dr r 
with solution 7 
B= 2.14 
r+C ( ) 


where C is a constant. Inserting this in (2.10) then yields a first order equation for A 
of which the solution is 
a(r+C) 


r 


A= (2.15) 


where a is a further constant. This latter constant effectively just changes the units 
of time, and it is sensible to choose this so that the speed of light c is 1, which we 
temporarily employ. We have thereby recovered the Schwarzschild metric, (2.7), and 
can identify the constant C as —2G M. 

One notices straightaway that some kind of singularity exists at r = —C = 2G'M, 
but it is worth noting that this is not a singularity of the Riemann tensor (the only 
non-zero quantity transforming as a full tensor we have available for investigation). 
The entries of this all behave with r like 1/r?, corresponding to tidal forces which 
make neighbouring particles move apart or together in their motion, and none have a 
singularity at r = 2G'M. Indeed, the form in which Schwarzschild first discovered his 
metric also had no singularity at this radius (see e.g. and the comments in ÎTI), 
and this coupled with lack of singularity of the Riemann tensor except at the origin, 
perhaps contributed to the uncertainty stretching over many years as to the physical 
status of the distance 2GM. 

Nowadays, we of course recognise this as the “horizon”, the point where light be- 
comes trapped. We discuss this in more detail below, particularly in connection with 
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the Kerr solution, but we can understand this in simple terms by looking at the ‘coor- 
dinate speed’ of a light pulse. For such a pulse, the interval ds = 0, and for outward 
radial motion this means 
2 2 . dr A 
Adt” — Bdr =0, i.e. u` VB (2.16) 

The ‘horizon’ is therefore where the metric coefficient B goes to infinity, since this 
marks the point where light is no longer able to move outwards. This happens in the 
Schwarzschild case at r = rs =2GM/ 2, which is the location of the event horizon. 

Before proceeding further, we mention a few physical facts about Schwarzschild 
black holes. The first is that the radius of the event horizon, rg, corresponds to 3km 
per Solar mass. This means that the mean density of a black hole is given by 


M —2 
p=2x 10! GE) gmce™!. (2.17) 
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Thus black holes with masses above 108 Mọ, have average densities below that of 
water, or the Sun. Those above a few billion Mo have densities below that of air. So 
from a mean density point of view, supermassive black holes are not particularly dense. 
The light-crossing time of the event horizon (i.e. a length equivalent to its diameter) is 
0.2 ms for a 10 Mọ black hole. 20 s for 10° Mo and about one day for 5 x 10°? Mọ. 


3 Relativistic Astrophysics emerges 


The previous neglect of black holes changed in 1963 when two important events oc- 
curred, the discovery of quasars by Maarten Schmidt and the discovery of the 
solution for a spinning black hole by Roy Kerr [BT]. These combined to lead to Rela- 
tivistic Astrophysics as the combination of Astrophysics and GR, embodied in the first 
Texas Symposium of Relativistic Astrophysics in December 1963. Kerr’s brief talk 
has in retrospect been called the most important announcement at the Symposium but 
was not mentioned by any of the three summarizers of the meeting [62]. 

The discovery of quasars resulted from the accurate position for the cosmic radio 
source 3C273 obtained by Cyril Hazard and collaborators using the lunar occul- 
tation technique. Schmidt used the 200" Hale Telescope to take an optical spectrum of 
the starlike object at that position in the Sky. Some weeks later he identified emission 
lines in the spectrum with Balmer lines of hydrogen, redshifted by 0.16. Interpreting 
the redshift as due to cosmic expansion means that it is at a (luminosity) distance of 
0.75 Gpc. From its apparent brightness, the object can be inferred to be more than 
10!? times more luminous than the Sun. 3C273 can be found on optical plates taken 
since the late 1800s and was soon shown to be variable, at times changing significantly 
within a week. Ignoring relativistic considerations, causality requires that significant 
variations cannot occur on timescales shorter than the light crossing time of the object. 
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It therefore appears that the prodigious power output of 3C273 emerges from a region 
of size similar to the Solar System! |! 


Figure 1. The first quasar identified, 3C273, as imaged by the Hubble Space Telescope. 
It lies over two billion light years away and is the brightest quasar in the visible Sky. The 
structure to the SW (lower right) is the outer parts of its relativistic jet, which is seen 
from radio to X-ray wavelengths. The projected length of the jet is over 200 thousand 
light years, so its true length is nearly an order of magnitude larger than our galaxy. 
The mass of its black hole has been measured by the optical reverberation technique at 
almost a billion Solar masses. The X-shape is an instrumental artefact created by support 
structures in the telescope; it implies that the source is unresolved. 3C273 varies on 
week-long timescales indicating that the emission region is less than a light week across. 


Kerr’s solution to the Einstein field equations was not immediately recognised as 
the exact solution for a spinning black hole. We now know the geometry described is 
the unique and complete description of the external gravitational field of an uncharged 
stationary black hole. It is now thought that the power of 3C273 and other quasars is 
due to accretion onto a Kerr black hole, which we now give some preliminary details 


' 3C273 is not just a pointlike object but has an associated linear structure which we now know to be 
a relativistic jet, but the above argument is good for the main central source. 


Astrophysical Black Holes 7 


of. 


3.1 Rotating black holes 


It is now known that the totality of black holes can be described using just three pa- 
rameters. Whatever makes them up, and however they were formed, in the end only 
three parameters matter — their mass, M; charge, Q and angular momentum J. Due 
to the high conductivity of interstellar matter, black holes are not likely to have a net 
charge for long, so the only relevant parameter in addition to mass is spin] A metric 
applicable to a black hole with angular momentum was found by Kerr in 1962 [31]. 
This metric, in the form later developed by Boyer & Lindquist 6], is 


Fee. (1 | ip 1 [eresie 
c 


P 2 
t H 
5a Z dtdo Agr 


(3.1) 


Ten? ein2 
+ pd? 4 G ta? 4 eas = *) sin? sag], 
pc 


where a = (J/Mc) is the angular momentum of the black hole per unit mass (and 
has the dimensions of distance), A = r? — (2GMr/c*) + a? and p = r? + a? cos? 0. 
Substituting this metric into equation for the connection coefficients and calcu- 
lating the Riemann from these, the Ricci tensor Ryy = R> pv will be found to 
vanish, and therefore satisfy the vacuum field equations (2.3). This is a formidable 
calculation, however, and techniques which derive the metric in the context of more 
general azimuthally symmetric spacetimes (e.g. ) are in fact easier to carry out. 

Note if the black hole is non-rotating, then J = a = 0 and the Kerr metric reduces 
to the standard Schwarzschild metric (2.7). Almost certainly, all real black holes in 
the universe will be of the Kerr type. The idea is that just as infalling matter will have 
angular momentum, so will the material from which the black hole formed, leaving it 
both with a mass and net angular momentum. 


3.2 Black holes as energy sources 


The enormous luminosity of 3C273 and other quasars was soon linked to black holes 
by Edwin Salpeter [60] and (separately) Yakov Zeldovich [78]. The energy is due to 
the deep gravitational potential well of the black hole which can be liberated by colli- 
sions outside the event horizon as gas falls in. If two gas clouds with small amounts of 
equal but opposite angular momenta fall into a Schwarzschild black hole colliding at 
4r, (note frg is defined as GM/ cœ) and momentarily coming to rest with their kinetic 
energy emitted as radiation, then about 29% of the rest mass energy of the clouds is 


> Newman et al found the solutions to the Einstein-Maxwell field equations in 1965 and the 
resulting Kerr-Newman geometry describes spinning charged black holes. 
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released. More generally the angular momentum will be much larger and the radiative 
efficiency of accretion smaller. A specific model involving an accretion disc was pro- 
posed by Donald Lynden-Bell in 1969. Due to the small size of a black hole, it 
is most unlikely that gas falling into a black hole does so radially, but will have some 
angular momentum which will cause it to orbit the hole. Viscosity in the swirling 
gas will cause matter at smaller radii, which is moving faster, to transfer its angular 
momentum outward to material at larger radii, which is moving slower. The gas then 
spreads in radius forming an accretion disc in which angular momentum is transferred 
outward as matter flows inward. The accretion disc model of Lynden-Bell was then 
studied in a detail in the early 1970s by Pringle & Rees [58], Shakura & Sunyaev 
and, for the Kerr metric, by Novikov & Thorne [49]. The gravitational energy lib- 
erated by the inflow heats the disc which radiates locally as a quasi-blackbody. The 
disc is thin but may extend outward for a considerable distance. This basic picture 
probably accounts for much of the energy liberated by accreting black holes. There 
are important modifications due to the magnetic nature of the ionized infalling plasma 
which will be discussed later. Accretion onto a black hole is the most mass-to-energy 
efficient process known, apart from direct matter-antimatter annihilation which, due 
to the rarity of antimatter in our universe, is highly uncommon. Such accretion may 
account for 20-30% of the energy released in the Universe since the recombination era. 
To understand the details of how this can happen, we now examine particle motion in 
GR, first generally, and then applied specifically to the question of efficiency of grav- 
itational energy release around a Schwarzschild black hole. (The more complicated 
issue of energy release around a Kerr black hole is treated in Section|7.1]}) 


3.3 Motion in the Schwarzschild metric 


The key to most astrophysical applications of the Schwarzschild metric is how point 
particles and photons move in it. In General Relativity, particles move on geodesics of 
the metric, i.e. the paths with an extremal lapse of proper time (for a massive particle) 
or ‘affine parameter’ (for a massless particle), along the worldline. If we let ds be the 
differential interval along a path, and s 4p the total interval between two points A and 
B ona given path, then with &,, = dx,,/ds we have 


a B 1 a dx” dx” 
SAB L ds I [guu dz" dx” |? I ow ds =| ds 


B 
=l G(x”, ś")ds where G(x", 4") = [gurit i”]2 
A 


Nie 


(3.2) 


and finding the path which extremises s 4p then leads to the Euler-Lagrange equations 


(one for each u) 
d {0G OG 
ds (a) B ee 


Astrophysical Black Holes 9 


Note using the result is is easy to demonstrate the equivalence of with the 
alternative, and perhaps more familiar form 

da | Te dx? 0 (3.4) 
ds? ds ds 
The advantage of the former, B3}. comes from the fact that it does not need knowledge 
of the connection coefficients to compute it, and also conservation laws for variables 
of which the metric is not an explicit function, are easy to read off. Note from its 
definition (compare the second and fifth quantities in the chain of equalities ), 
that as a numerical quantity G(x”, <") evaluates to 1 (at least for a massive particle), 
which we can use to simplify formulae once its functional dependence has already 
been used. 

So we now apply these results to the metric in the general static form given by 
2.7). (The advantage of carrying this out for the general form, rather than just 
Schwarzschild, is that it enables us to consider results for other types of black hole, 
such as Reissner-Nordstrom, and Schwarzschild-de Sitter — see below.) We can as- 
sume w.l.o.g. that the particle motion is confined to the 0 = 1/2 plane, and then the 
equations we find are 


AË — Bi? -r° =1 (3.5) 


coming from G(x”, t”) = 1, and 
Ai=k, and rd=h (3.6) 


coming from the Euler-Lagrange equations in t and ¢ respectively, and where k and h 
are constants. These last two equations result from the constancy of the metric coeffi- 
cients in t and ¢ and correspond to the conservation of energy and angular momentum. 
For a particle of mass m, kmc* can be identified as the particle energy, and h is the 
specific angular momentum per unit mass. 

There is no point forming the Euler-Lagrange equation in r, since we already have 
a simple expression for ù available from the combination of and (3.6), yielding 


LIke h 
r=5(5-5-10) (3.7) 


We return to this general form later, but now wish to obtain results specific to the 
Schwarzschild case. With A = B! = (1 — 2GM/r), and reinstating c for the 
remainder of this section, we obtain 


2 


H : 3.8) 

r2 ' ( 
Using the usual Newtonian substitution u = 1/r, and changing the independent 
variable to azimuthal angle ¢ rather than the interval s, so that it is easier to discuss 
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the shape of the orbit, we differentiate to obtain 


u GM 3GM , 
ip u= 72 | z u. (3.9) 
In Newtonian gravity, the equation for planetary orbits, in the same notation, is 
Qu GM 


so we have almost got the Newtonian answer, except for the extra term 3GMu? / e. 
This is what gives all the relativistic effects, and we note it correctly goes to zero as 
c> œ. 


3.4 Circular orbits 


An important application of the orbit formula in the context of high energy astro- 
physics, is what it tells us about circular orbits in Schwarzschild geometry. These will 
be approximately the orbits of material accreting onto black holes, since infalling ma- 
terial nearly always has angular momentum, and we would not generally expect direct 
radial infall. 

If r is constant, then our equation for u yields 


Mr? 
h? = GME: (3.11) 
Putting è = 0 in equation gives us 
aet a (3.12) 
Cr pe nee 


and then putting both these last two results together yields an equation for k in terms 
ofr alone. We derive: 


i 2GM 
2 
A (3.13) 
_ 3GM_ 
re 
Remembering k = E/m@, where E is the particle energy, we have found that the 
energy of a particle in a circular orbit is 
_ 2GM 
Pg 
Ecirc = = me" (3.14) 
_ 3GM 
re 
An obvious check on this equation, is whether it can reproduce the Newtonian expres- 


sion for the total energy of a circular orbit in the limit of large r. Using the binomial 
theorem we see that indeed the first two terms in an asymptotic expansion in r are 


GM 
Bar ~me — i aeS 
2r 


(3.15) 
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Thus we get agreement at large r with the usual Newtonian expression (derived via 


Í M -GM 2 M 
By =e ee = mm GMm = GMm if mu G mi (3.16) 
2 r 2r r r2 


provided we realise that it enters as a correction to the rest mass energy mc’, which is 
the dominant term. 

The equation we have just found for the energy of a circular orbit, provides us with 
a remarkable amount of information about the nature of such orbits. 

First we see that in the limit m — 0, the orbit r > 3GM/c is of interest, since 
the singularity in the denominator can cancel the zero at the top. In fact this is the 
circular photon orbit at r = 3GM/c’, which we will derive below in a treatment of 
photon motion. 

Secondly, we can see which orbits (for particles of non-zero rest mass) are bound. 
This will occur if Eare < mc?, since then we have less energy than the value for a 
stationary particle at inifinity. The condition for Eire = mc? is that 


=1- -5 (3.17) 


(: E a _, 3GM 


re 


This happens for r = 4GM/e orr = œ. Thus over the range 4 < r < o the 
circular orbits are bound. 

This appears to show we can get as close as 2 Schwarzschild radii to a black hole 
for particles in a circular orbit, suggesting that this is where the inner edge of any 
accretion disc would be. But is such an orbit stable? We discuss this first in the specific 
context of the Schwarzschild metric, and then later look at stability in the context of 
more general metrics of the form ; 


3.5 Stability of circular orbits 


In Newtonian dynamics the equation of motion of a particle in a central potential is 


2 
(5) +V(r) =E, (3.18) 


where V (r) is an “effective potential”. For an orbit around a point mass, the effective 
potential is 

h? GM 
Qr2 r? 
where h is the specific angular momentum of the particle. Since the 1/r? term will 
eventually always exceed the 1/r term as r — 0, we can see that in Newtonian dynam- 
ics a non-zero angular momentum provides an angular momentum barrier preventing a 
particle reaching r = 0 — see Fig.|2| In this effective potential, bound orbits have two 
turning points and a circular orbit corresponds to the special case where the particle 


V(r) = (3.19) 


12 A.C. Fabian and A. N. Lasenby 


=«——— effective potential 


unbound orbit 


elliptical orbit 


circular orbit 


Figure 2. The Newtonian effective potential showing how an angular momentum barrier 
prevents particles reaching r = 0. 


sits at the minimum of the effective potential. However, as we have already partially 
seen in equation (5.2). the same is not true in General Relativity. 
Starting with equation we can rewrite this as 


1, k 2GM\ GM 1, 
=-7C 


2 Op? Cr r 2 


(k? — 1) (3.20) 


where we recall k = Ean /mc* andr? =h. 
Thus, although the r.h.s. is not the particle energy here, the fact that it is constant 
tells us that 


U(r) (3.21) 


— 2r? 7 
is an “effective potential” for the problem, which we can use to study stability in the 
same way as in the Newtonian case. Note that the relativistic term (1 — 2GM/c’r) 
weakens the centrifugal effect of angular momentum at small r. 

Differentiating this expression, 


dU i? 3VGM _ GM 


CT E 


B = ( ge) GM 


i (3.22) 


dr r3 


and so the extrema of the effective potential are located at the solutions of the quadratic 


ert r2 


equation 
3RGM 


GMr* —h?r4 5 
c 


0, (3.23) 
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i.e. at 


h? GM\* 


Ifh=v 12¢¢ there is only one extremum, and there are no turning points in the orbit 
for lower values of h. At this point r = 6GM/c? = 3R,. Fig.|3| shows the effective 


Figure 3. The effective potential U(r) plotted for several values of the angular momen- 
tum parameter h (units here have GM/c? = 1). 


potential for several values of h. The dots show the locations of stable circular orbits. 
The maxima in the potential are the locations of unstable circular orbits. 
What is the physical significance of this result? The smallest stable circular orbit 
= 6GM 
Tmin = — 7: (3.25) 


Gas in an accretion disc settles into circular orbits around the compact object. How- 
ever, the gas slowly loses angular momentum because of turbulent viscosity (the tur- 
bulence is thought to be generated by magnetohydrodynamic instabilities). As the gas 
loses angular momentum it moves slowly towards the black hole, gaining gravitational 
potential energy and heating up. Eventually it loses enough angular momentum that 
it can no longer follow a stable circular orbit and so it falls into the black hole. On 
this basis, we can estimate the efficiency of energy radiation in an accretion disc via 
looking at a plot of the ‘fractional binding energy’ E/(mc*) —1 = k — 1 versus r (see 


Fig. H). 
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Bound energy for circular Schwarzschild orbit 


(E/(me) -1)% 


r/p 


Figure 4. A plot of E/(mc*) — 1 in percent versus r (the latter measured in units of 
u = GM/c’) where E is the energy of a particle of mass m in a circular orbit at radius 
r about a Schwarzschild black hole. 


The maximum efficiency is of order the gravitational binding energy at the smallest 
stable circular orbit divided by the rest mass energy of the gas. From the plot we can 
see that this will be about 6%. More accurately, from equation we see that at 
r=6GM/c k = E/(mc’) is 2\/2/3, hence we obtain for this efficiency 


Eace © 1 — 22/3 ~ 5.7% (3.26) 


The equivalent Newtonian value, is not far away at 


1GMm 1 1 
N ~ ~ 2 27 
cal Tmin Mce 12 A ere) 


As we will see, this value can be even larger for a black hole with spin, and an accre- 
tion disc can convert 5-20 percent of the rest mass energy of the gas into radiation, 
depending on spin. This can be compared with the efficiency of nuclear burning of 
hydrogen to helium (26 MeV per He nucleus), 


€nuclear ~ 0.7% (3.28) 


Accretion discs are capable of converting rest mass energy into radiation with an effi- 
ciency that is about 10 times greater than the efficiency of nuclear burning of hydrogen. 
The ‘accretion power’ of black holes causes the most energetic phenomena known in 
the Universe. 
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4 Evidence from X-rays, quasars and AGN 


The early 1960s also saw the birth of X-ray astronomy, when in June 1962 Riccardo 
Giacconi and colleagues launched a sounding rocket with Geiger counters sensi- 
tive to X-rays. As the detectors scanned across the Sky during their 10 minutes above 
the Earth’s atmosphere, a broad peak of emission was discovered from the general di- 
rection of the Scorpius constellation and a steady background of cosmic X-rays were 
seen. Sco X-1 and the cosmic X-ray Background had been discovered. The rapidly 
variable X-ray source Sco X-1 is a neutron star accreting matter from an orbiting com- 
panion star and the X-ray Background is now known to be the integrated emission 
from accreting supermassive black holes at the centres of distant galaxies (see Fig.|5). 


Figure 5. The XMM COSMOS field is about 30 arcmin across and shows the X-ray 
Background resolved into pointlike and extended X-ray sources. The pointlike sources 
are mostly accreting black holes (distant AGN) and the extended sources are due to hot 
gas pooled in the gravitational potential wells of clusters and groups of galaxies. 


The first clear association of X-ray astronomy and black holes was made after Giac- 
coni and colleagues studied the brightest X-ray source in the constellation of Cygnus, 
Cyg X-1, using the first X-ray astronomy satellite, Uhuru, which continuously scanned 
the Sky. Cyg X-1 showed chaotic variability on timescales down to a fraction of a sec- 
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ond. It has two relatively long-lived states, one when the spectrum is soft and the other 
hard. A radio source was found to appear when it switched into the hard state and 
accurate measurement of the position of that source enabled it to be identified with a 
7th magnitude B star, HDE 226868. Separate teams of optical astronomers, Webster & 
Murdin and Bolton [5], then found that this massive star was being swung around 
at 70kms~! by an unseen object on a 5.6 day period. This enabled mass estimates to 
be made (which relied on understanding the mass of the companion B star) yielding 
at least 3.5 which is above the likely upper mass limit for a neutron star of about 
3Mo. The small size implied by the rapid variability combined with the lower mass 
limit pointed to a black hole in Cyg X-1. 

Many other black hole X-ray binary systems are now known, some of which have 
much clearer mass estimates. M33 X-1, for example, in the nearby galaxy M33 has 
a precise distance and stellar mass estimate and also shows eclipses which enable the 
orbital inclination to be deduced, leading to a mass for the X-ray source of 15.64 + 
1.45Mo BI]. 

By the 1980s, many quasars and X-ray binary systems were known, showing that 
black holes are common. Quasars were seen to be an extreme part of the more general 
phenomenon of Active Galactic Nuclei (AGN). The centres of a few percent of all 
galaxies appear to have some non-stellar activity in the form of bright broad emission 
lines, a nonthermal radio source or an X-ray source. Quasars occur when the AGN out- 
shines the host galaxy in the optical band. AGN of lesser power are known as Seyfert 
galaxies or just low luminosity AGN as the power drops. Jets of highly-collimated rel- 
ativistically outflowing plasma are found in about 10% of AGN. An example outflow 
in one of the nearest AGN is shown in Fig. (6| 


Unusual emission spectra at the centres of some galaxies had been known as a phe- 
nomenon since about 1908, and the first jet was reported in M87 by Heber Curtis in 
1918 [T2]. Little follow-up work was done before the 1960s, apart from Carl Seyfert’s 
PhD thesis in the 1940s [64]. By the end of the 1980s most astronomers considered 
that AGN were powered by accretion onto black holes, but a minority still favoured 
some other explanation, such as multiple supernova outbursts. The arguments for black 
holes were strong but circumstantial. A major observational problem was that black 
holes are by definition difficult to observe directly. 

Quasars appear to be associated with a past phase in the evolution of the Universe, 
when it was 1-5 billion years old (mostly at redshifts of 2-4). They are now relatively 
uncommon compared with that era. (3C273 is a rare quasar at low redshift: low is rela- 
tive here as it was moderately high when discovered in the early 1960s.) The enormous 
powers involved mean that the accretion rates were high and the mass doubling time 
of the black holes could have been a few 100 million years. Radiation pressure would 
have restricted the inflow to below a limit first deduced for stars by Eddington, known 
as the Eddington limit. The limit is obtained by comparing the force of radiation acting 
on an electron (through Thomson scattering of cross-section or) at the surface of the 
object, with the force of gravity on a proton there. The electron and proton, although 
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Figure 6. Centaurus A is one the nearest radio-loud AGN. The image on the left is a 
combination of infrared through X-ray images, its jets are seen going from the NE to 
SW. On the right is a larger scale image of low frequency radio emission showing the 
large diffuse lobes of radio-emitting plasma ejected by the accreting black hole out to 
intergalactic space. 


assumed free, are bound together by electrostatic attraction. 


GMm, L hv 

R TIT aap e” 
where the radiative force term is the flux of photons (of typical energy hv) times the 
cross-section times the momentum of a photon. This gives 


A4rGM m,e 
Lea = ————>,, (4.2) 
OT 


(4.1) 


F, grav — 


which is about 10°! W per Solar mass. 

Since the Eddington limit is proportional to the mass of the black hole the black 
hole mass can grow exponentially, if sufficient fuel is available. If growth is at the 
Eddington limit, the e-folding time of the black hole mass is about 4 x 107 yr. This 
timescale is sufficient for the typical massive black hole to have grown from stellar 
mass (say 30M.) but it does become challenging for the most distant billion Solar 
mass quasars found above redshift 6. 


5 Sgr A*, the Black Hole at the Centre of the Milky Way 


The evidence for the existence black holes changed though the 1990s owing to careful 
observations of the motion of stars moving around the centre of our own Milky Way 
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galaxy, first starting in 1991 by Reinhard Genzel and colleagues using ESO tele- 
scopes and later joined by Andre Ghez and colleagues using the Keck telescopes. 
An interesting picture of the technique used to achieve the required image stability for 
this work is shown in Fig. [7] 


Figure 7. The central stellar bulge of our Galaxy. A laser beam from one of the ESO 
Very Large Telescopes points at the very centre where Sgr A* resides. This creates an 
artificial star in the Earth’s upper atmosphere which is used to stabilise images made of 
the stars orbiting the central black hole. 


Most stars, like the Sun, orbit the centre of the Galaxy at about 230kms~!, but 
within the innermost light year they move faster until at a distance of 2000r, from the 
dynamical centre a star orbits at up to 3000 kms~!. The inferred mass within that orbit 
is 4 million Mo and the dynamical centre coincides with a strange radio source in the 
constellation of Sagitarius long known as Sgr A¥*. It is also a flickering infrared and 
X-ray source. No stellar emission which can be attributed to say a massive cluster of 
normal stars is seen. The stellar orbits of dozens of stars there are known in 3D, both 
from their apparent motions mapped in the plane of the Sky and from line-of-sight 
Doppler shifts. The orbit of one of the closest stars, S2, is shown in Fig. [8] [26]. The 
mass density within the innermost orbits exceeds 10!8 Mo pe~? and there is nothing 
known to Physics which can lie there other than a black hole. Clusters of stars, even 
neutron stars, would interact and collide. The above lower limit on the mass density 
of the central object is increased by a further 5 orders of magnitude when constraints 
on the maximum size of Sgr A* obtained from millimetre band Very Long Baseline 
Interferometry (VLBI) are included. 

The origin of the young, massive stars which are bright enough to be tracked around 
the Galactic Centre is not yet fully understood. The underlying stellar density there is 
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The orbit of S2 (1992-2013) 
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Figure 8. The path of star S2 about SgrA*. The elliptical orbit has a period of 15.2 
years, a major axis of 5.5 light days and is inclined at 46 degree to the plane of the Sky. 
It requires that a mass of 4 x 10°Mo lies within a radius of 17 light hours. The only 
long-lived object which is physically consistent is a massive black hole. 


so high that it should inhibit further star formation. 

As well as the issues of mass within a given radius, there is also hope that general 
relativistic effects on the orbits of stars and possibly gas clouds near Sgr A* will be 
observed in the near future, so we discuss now the details of these effects, including 
the important aspect of ‘capture’ by a black hole. 


5.1 GR effects on orbits 


We derived above the ‘shape’ equation for orbits around a Schwarzschild black hole, 
3.9). For the Newtonian case, then due to the harmonic kernel at the left, and the 
constant at the right, it is obvious that solutions of the form u = a+ bcos ¢ will work, 
i.e. ellipses, so if the extra term in the relativistic version is small we expect the 
orbits to be modified ellipses. In the Solar System, in which a Schwarzschild metric 
due to the Sun applies, these modifications are very small (e.g. the largest effect is 
the 42 seconds per century precession of Mecury’s orbit), but in the Schwarzschild 
metric around a black hole, we can expect much larger effects. For example, in Fig. 9 
we show the motion for an object starting about 10 Schwarzschild radii out in an 
almost circular orbit. The gaps between the markers on the orbit are laid down at 
equal intervals of time, and so indicate how fast the object is moving. We can see that 


20 A.C. Fabian and A. N. Lasenby 


strongly precessing ellipses are obtained, but with a motion that looks as though it will 
continue indefinitely (which it does). 


However, in Fig. we show what happens when starting from the same point, and 
also moving tangentially, but now with a smaller specific angular momentum. This 
time, the object is clearly ‘captured’ by the black hole, and ends up crossing the black 
hole horizon. 
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Figure 9. Advance of perihelion in a Schwarzschild metric. The units of distance are 
GM_/c*, where M is the mass of the central body. The particle is launched tangentially 


and given a specific angular momentum h of 3.75 where the circular orbit h would be 
4.85. The circle marks the Schwarzschild radius of the central black hole. 


We can gain some insight into how the ‘capturing’ comes about, by considering the 
orbit equation in a different form. For a massive particle, the interval s we have been 


using is just the proper time of the particle, which we label T. If we start with equation 
(3.8), differentiate with respect to interval T and then remove first derivatives, dr /dt 


using the original equation again, one arrives at the following: 


@r GM R 


r2 r3 


3h°GM 
<a (5.1) 
Cr 
The first two terms are very Newtonian-like, and correspond to an inward gravita- 


tional force and a repulsive term, proportional to angular momentum squared, which 
is basically the ‘centrifugal force’. What is new is the third term, also proportional 


to angular momentum squared, but this time acting inwards. This shows that close to 
the hole, specifically within the radius r = 3G.M/c, centrifugal force ‘changes sign’, 


and is directed inwards, thus hastening the demise of any particle that strays too close 
to the hole, as in the example of Fig. 
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Figure 10. As for preceding figure but now for a particle projected with h = 3.5 (in units 
with GM/c? = 1). 


Orbital precession 


The equation we wish to solve is 


@u GM ,3GM » 


dd? ru k` g (u 


1 
= -) (5.2) 
r 


in the limit that the departure from Newtonian motion is small. This would apply to 
the motion of the planets in our solar system for example. 
The Newtonian solution to this equation is 


u= (M iero (5.3) 


h2 
and so we can use this as a first approximation, and then iterate to get a better one. 
Substituting into the r.h.s. of (5.2), we obtain the new equation 


du GM 3(GM) 


{e SG ee 


(1 + ecos ġ}. (5.4) 


This will be solved by u = Ge + ecos ġ) + the particular integral (P.1.) of the 
equation: 

Pu 79 

ag a ener cos” ġ) (5.5) 
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where A = 3(GM)3/c7h’ is very small. The PI. can be found to be: 


A(I+ebsing +e (= - ze0s20) ) : (5.6) 


Now, in this expression, the first and third terms are tiny, since A is. However, the 
second term, Aeg sin ¢, might be tiny at first, but will gradually grow with time, since 
the @ part (without a cos or sin enclosing it) means it is cumulative. We must therefore 
retain it, and our second approximation is 


GM 
u= yC + ecos ġ + dedsin ġ) (5.7) 
where ; 
3(GM 
ô= our <1. (5.8) 
Using 


cos (¢(1 — ô)) = cos ġ cos ôo + sin ġ sin ô ~ cos o +sing ôo for ô<« 1, (5.9) 


we can therefore write 


GM 
ux pz (+ ecoslol — ô)]). (5.10) 
u is therefore periodic, but with period Ka The r values thus repeat on a cycle which 
is slightly larger than 27, and we find 


2 
ze 2r 276 = OM 


Ad = I5 M (5.11) 


But from the geometry of the ellipse, and the Newtonian solution, where we know 
l= h? /GM and l = a(1 — e”), we can get the final result: 
_ 6n(GM? — 6nGM 

eI(GM) a(l = e) 


Ad (5.12) 


For example, Mercury’s orbit has a = 5.8 x 10!° m, eccentricity e = 0.2 and we know 
Mo =2 x 10°°kg. Therefore our prediction for the precession is 


A@ = 5x10” radians per orbit 
0”.1 per orbit. 


Since Mercury’s orbital period is 88 days, we would thus expect to accumulate a pre- 
cession of 43" per century. This is what is observed after correction for the perturba- 
tions due to the other planets, which cause a total precession of more like 5000" per 
century. 
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Turning now to the Galactic Centre, the star S2, whose orbit was shown in Fig. 
has an eccentricity of 0.876 and semi-major axis of 980 AU [I6\. Taking the mass of 
the Galactic Centre black hole as 4 x 106 Mọ, then our formula yields a precession of 
11 arcmin per orbit, which takes about 15 years. This sounds large, but projected on 
the sky at perihelion amounts to only about 0.5 milliarcsec. By the use of near-infra red 
interferometry, this may be observable in the relatively near future (see Section (70), 
but as for the solar system it is competing effects from other nearby masses, and the 
general matter distribution near the centre, that will be key to determining whether it 
is the GR effect itself which is being seen [44]. 


6 AGN Feedback 


It was earlier unclear whether galaxies with central black holes were special in some 
way or not. Quasars could be long lived in just a few galaxies or short lived and occur 
in every galaxy. Then careful imaging of the centres of nearby galaxies with HST and 
other telescopes revealed through the 1990s that most galaxies host a massive black 
hole. The mass of the black hole is correlated with the mass of the galaxy bulgd’|that 
it is embedded in (the Magorrian relation [37]). Current results indicate that the host 
bulge has a mass 500 times that of the black hole. Some argue that the correlation is 
better if the velocity dispersion of the surroundings stars (beyond the radius where the 
black hole’s mass dominates) is used (the gy — ø relation — an example of this is 
shown in Fig. [IT). 

Both correlations show considerable scatter but hold over several orders of magni- 
tude. 

The M —o relation poses the question of how the black hole at the centre of a galaxy 
knows the total mass of the host galaxy? The ratio of the physical size of the black hole 
event horizon to the galaxy can be 100 million or more (a similar scaled comparison is 
the Earth and a football). Some astronomers argue that larger galaxies just have larger 
black holes, but the correlations look better then that. An exciting possibility is that 
the question should be turned round and it is the black hole which controls the mass of 
the host galaxy. In 1998 Silk & Rees proposed that it was the phenomenal energy 
output in the quasar phase, as the black hole grew in mass by accretion, which was 
responsible. Both grow from gas in the galaxy — the galaxy forming stars from the 
gas, the black hole accreting the gas. Eventually the black hole is so big and powerful 
that it blasts all the gas from the galaxy so both stop growing and we are left with the 
observed relation. 

There is clear evidence that massive black holes have grown by accretion from the 
Soltan relation [67]: 


Excell +2) = T p’ (6.1) 


> This means that the mass of any galactic disc is ignored. 
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Figure 11. The black holes mass — stellar velocity dispersion relation Myn — o (27). 


where z is the mean redshift at which the accretion occurs. It is a recasting of the 
famous equation, E = mc?, in terms of densities. The factor of 1 + z is due to 
the redshifting of the energy of the radiation; there is no such factor for mass. 77 is 
the radiative efficiency of the accretion process L = nMc*. The energy density of 
radiation from accreting black holes, Eacc can be measured from the summed spectra 
of AGN (and the X-ray Background) and pgy can be estimated from the mass function 
of galaxies together with the Mga — Mg, relation. Results show agreement if n is 
about 0.1 which, as discussed later, is typical for luminous accretion. 


It is straightforward to show that the growth of the central black hole by accretion 
can have a profound effect on its host galaxy. If the velocity dispersion of the galaxy is 
a then the binding energy of the galaxy bulge, mass Mga, is Egal © Mgao’. The mass 
of the black hole is on average observed to be Mpy ~% 2 x 1073 Mgal: If the radiative 
efficiency of the accretion process of 10%, then the energy released by the growth of 
the black hole is given by Egy = 0.1 Mpu. Therefore Epy/Ega © 2 x 1074 (c/o}?. 
Very few galaxies have ø < 350kms™!, so Egy/Ega > 100. The energy produced 
by the growth of the black hole therefore exceeds the binding energy of the host galaxy 


Astrophysical Black Holes 25 


by about two orders of magnitude! 

If even a small fraction of the accretion energy can be transferred to the gas, then an 
AGN can have a profound effect on the evolution of its host galaxy. In practice, radi- 
ation from the accreting black hole will not have any significant influence on existing 
stars in the galaxy, It can however strongly influence gas clouds from which new stars 
can potentially form. By ejecting those clouds from the galaxy, AGN feedback can 
effectively stop further stellar growth of a galaxy. 

The original mechanism of Silk & Rees is actually a little too effective and although 
it predicts a slope to the relation which is acceptable, the normalization is too low. Gas 
is blasted out too easily. However, energy is likely to be radiated away in the process 
which is likely to be controlled by conservation of momentum 45]. (Rocket 
science would be much easier if only energy conservation were at stake.) It has since 
been shown that this leads to the correct normalization and a slightly flatter slope. The 
details now centre on the precise mechanisms responsible for the ejection (e.g. 
radiation pressure or accretion disc winds). Energetic winds are seen from some AGN 
and UV radiation interacts strongly with dusty gas in all quasars exerting an outward 
force through radiation pressure. 

It is interesting to note that the effective Eddington limit for dusty gas exposed to the 
ultraviolet radiation of a quasar is about 500 times that of ionized gas such as expected 
close to the quasar. The M — ø relation means that when the quasar is at its Eddington 
limit locally, the host galaxy is also at its effective Eddington limit globally. 

Dramatic evidence of AGN feedback occurs in the cores of many clusters of galax- 
ies. The intergalactic medium in clusters has been squeezed and heated by gravity 
to create an intracluster medium which is 10 million K or hotter. The baryons in the 
hot gas exceed all baryons in the galaxy members of the cluster by a factor of more 
than ten. The most massive galaxies known lie at the dense centres of clusters where 
accretion onto their massive black holes feeds energy into the intracluster medium, 
preventing it from cooling down due to the emission of X-rays. The total mass of gas 
in a cluster and the depth of the gravitational well is so large that the gas cannot be 
expelled. Energy instead flows from the black hole accretion process into the gas to 
maintain a themodynamic status quo. 

The energy flows from the near vicinity of the black hole in the form of relativistic 
jets which push the hot gas back, forming giant galaxy-sized bubbles of relativistic 
plasma (cosmic rays and magnetic field) — an example is shown in Fig. [12 

The bubbles are buoyant and when large detach and float upward in the cluster 
potential, while new bubbles form and grow. The power in the bubbling process can 
be estimated from the energy content of a bubble (4PV where P is the pressure of 
the surrounding gas and V is the volume of the bubble, the factor 4 being appropriate 
for relativistic plasma) divided by the bubble growth time. This power matches the 
radiative energy loss of the inner cluster core, indicating that close feedback has been 
maintained for billions of years [11]. 

AGN feedback implies that black holes play a central role in the final growth and 


26 A.C. Fabian and A. N. Lasenby 


Figure 12. X-ray image of the core of the Perseus cluster centred on the giant galaxy 
NGC1275. The image is about 100 kpc or 300 thousand light years across. The AGN at 
the centre has blown 10 kpc diameter bubbles of relativistic plasma into the surrounding 
cluster gas. Detached buoyant bubbles are seen to the NW and S. Concentric ripples are 
seen in the hot gas due to the bubbling activity created by the AGN. It is clear that a black 
hole can have dramatic effects on gas structures both in and beyond the host galaxy. 


evolution of all massive galaxies. If black holes did not exist then massive galaxies 
would be much larger and brighter than we now see. 


6.1 Jets, Gamma-Ray Bursts and the Birth of Black Holes 


Relativistic jets occur in about 10% of AGN and in low-state BHB. They are thus a 
marker of the presence of a black hole. The details of how the material is accelerated 
and collimated are still unclear, with possibilities discussed later in this Chapter. Ap- 
parent superluminal motion indicates that the bulk Lorentz factor in many jets I ~ 10. 
Polarization indicates that the radio and much other emission are due to synchrotron 
radiation. The jets therefore contain electrons and a component with positive charge is 
expected, unless the jets represent enormous currents. Whether the positive particles 
are positrons or protons is not known. If they are protons (or even heavier nuclei), then 
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the power in the jets can be huge, exceeding the radiative power in some objects. 

Jets are presumably accelerated magnetically by the accretion disc. Fields emerg- 
ing from the disc are continuously wound up around and along the central axis. The 
details are not yet fully understood. Jets are seen from all classes of object exhibiting 
accretion discs, including young stars and accreting white dwarfs. A long discussed 
possibility from quasar jets is that in some cases the spin energy of the black hole 
is being tapped by the Blandford-Znajek mechanism, which is a variant of the Pen- 
rose process (see Sections [12.5] and [I2.7| below). There has been recent progress in 
successfully modelling jets using numerical simulations of this mechanism. 

The most relativistic jets occur in Gamma-Ray Bursts (GRB). These were discov- 
ered in the late 1960s by the Vela spacecraft which had a gamma-ray detector to mon- 
itor the Cold War Test Ban Treaty (but not made public by Ray Klebesadel and col- 
leagues from Los Alamos until 1973 (33). They are intense flashes of gamma-rays 
which have a flux which would make them readily visible to the naked eye if in the 
optical band. Many theories of their origin were proposed, mostly based on neutron 
stars. It was not until the late 1990s that a GRB was localized and identified with a 
distant galaxy, by the BeppoSax satellite, that made it clear that GRB are the most lu- 
minous events in the Universe, representing about a Solar rest mass of energy radiated 
within a few 10s of seconds. 

Luminosity and causality point to enormous concentrations of gamma-rays in a very 
small region which would instantly lead to electron-positron pair creation creating a 
electron scattering dense fog, preventing any radiation emerging. It is difficult to see 
how gamma-rays can be seen unless strong relativistic beaming, I ~ 100, is involved. 
Many aspects of GRB, including the radio through X-ray afterglow all support this 
highly relativistic interpretation [57]. The source of GRB must also be clean of any 
other matter, which could prevent the radiation emerging. It may then be surprising 
then that the source of GRB is now believed to be the very centre of dying massive 
stars as a black hole is formed. The swirling matter powers intense jets which burrow 
through and out of the star, to produce the GRB display. The birth of a black hole is 
therefore marked by one of the most amazing displays in the Universe. 

Stars which give birth to black holes are estimated to start life at masses exceeding 
25M. Such stars have luminosities which are a million times greater than the Sun 
and live for just a few million years. 


7 Current Observations of Accreting Black Holes 


Apart from the nearest examples, we know of the presence of most black holes from 
the power emitted as radiation and jets produced as a result of accretion. Matter falling 
into a black hole is accelerated by gravity up to the speed of light and if collisions 
occur then the matter can become very hot and radiate. The efficiency of the process 
depends on the rate of collisions and on whether the matter has time to radiate before 
falling into the event horizon. If the flow is a thin dense accretion disc then collisions 
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are plentiful and the efficiency depends on the ISCO (innermost stable circular orbit). 
Defining the radiative efficiency 7 through L = nMc’, then 7 ranges from 5.7% to 
42% as the spin a/M increases from 0 to 1. In practice it is impossible to spin a black 
hole up to exactly 1, and high spins mean 7 ~ 20 — 30%. Note that 7 for complete 
nuclear fusion (hydrogen to iron) is 0.7%, so black hole accretion is about an order 
of magnitude or more higher in efficiency. We now give some details of where these 
numbers come from, by looking at the orbits of massive particles around a Kerr black 
hole. 


7.1 Particle motion in the Kerr metric 


Before discussing the details of particle motion, we need to consider the general struc- 
ture of the Kerr solution, since this determines the character of the regions in which 
the particles are able to move. 


The most important regions for astrophysical purposes are shown in Fig. We 
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Figure 13. The stationary limit and outer horizon surfaces for the Kerr solution. 


have drawn an outer surface, called the ‘Stationary limit surface’, and the one lying 
inside this is called the ‘Outer event horizon’. Both of these are complemented by 
further surfaces (not drawn) lying closer to the black hole centre (the ‘Inner Stationary 
limit’ and ‘Inner horizon’ surfaces), but these are shielded by the outer horizon, and 
so do not have an immediate astrophysical role. 

The stationary limit surface marks the point where it is no longer possible for an 
observer to remain stationary in the (r,0,) coordinate system. Try as he or she 
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might, it is inevitable that they will be swept around by the rotation of the hole. Math- 
ematically, this corresponds to the point where the g Killing vector corresponding to 
invariance of the metric under time translations, changes from timelike to spacelike, 
i.e. to where the gu component of the metric passes through zero. This requires 
(setting c = G = 1 for clarity in the expressions) 


r? —2Mr+a’cos?@=0, ie r= M + vyv M? — a? cos? 0 (7.1) 


and the outer of these surfaces is shown. 

The other surface drawn corresponds to the point where the gry component of the 
metric becomes infinite, and marks the position of the event horizon. This is the 
position from beyond which a particle can never escape, whatever its mass or state 
of motion. It is different from the stationary limit, since particles in the inbetween 
region, although they cannot remain stationary, can in principle escape to infinity. 
This inbetween region is known as the ‘ergosphere’, since it is possible to arrange 
particle motions so that work can be extracted from it (see Section [72.5] below). The 
position of the horizon itself can be found by requiring 


A=r?—2Mr+a=0, ie r=M+VM?-o? io 


and again the outer of these two surfaces is shown. Given that this radius, unlike the 
stationary limit surface, does not depend on 0, it might be wondered why the surface 
is drawn oblate in the diagram. This is because that Boyer-Lindquist coordinates do 
not correspond to spherical polars even in the case where the mass vanishes, and 
an embedding of the 2d surfaces we are working with into Euclidean space (for the 
purposes of visualisation) results in ellipsoids even for r = const. — see Section 13.6 
of [BO] for details. 

The Lagrangian methods used above in the Schwarzschild case also work here, 
and one finds quickly the following results for the t and @ of a particle moving in 
the equatorial plane in terms of its conserved specific energy and specific angular 
momentum, k and h: 


k ((r +2M)a? + r°) — 2haM 
= rA 
h(r —2M) + 2kaM 
rA 


(7.3) 


d= 


As before, for the òè equation it is simpler to use the fact that evaluated numerically 
the interval function G(x", ¢”) is 1 for a massive particle, and 0 for a photon. We 
start with discussing a massive particle, since our initial emphasis is on the energy 
liberated from accretion disc orbits. Here we find 


2M 1 


t =k -14 (h? — a?(k* — 1)) 4 (h — ak)? (7.4) 
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We now have three equations we can integrate to find the motion. A key quantity for 
the energy and stability is the effective potential in r. We define this similarly to before 
via 


1 
a + Uker(r) = const. (7.5) 
resulting in 
M 1 M 
Ukerr(r) = T (h? — a?(k* — 1)) 3 (h — ak)’, (7.6) 


where an irrelevant overall constant has been ignored. 

If we compare with the equivalent in the Schwarzschild case, equation 3.27), we 
can see that despite the increased complexity of the Kerr metric, the effective potential 
has terms in just 1/r, 1/r? and 1/r? as before, resulting in overall similar behaviour 
(at least for small a). A difference, however, is that now the coefficients of the second 
and third terms depend on the particle energy as well as the angular momentum. 

For a circular orbit, we can use the fact that ù and ř are zero to find expressions for 
the specific energy and angular momentum. For the prograde orbits we are currently 
considering, these turn out to be 


1- tay% VMr — 72M +a 4 
k = , h= (1.7) 
y1- +2 4 y1- + 2a 5 


An important quantity for us to find is the radius of the innermost stable circular 
orbit, expected to be the inner edge of an accretion disc around a Kerr black hole, and 
the value the energy function takes there. 

As in the Schwarzschild case, we can use the effective potential to do this, asking 
that its second derivative is positive where the first derivative vanishes. Calculating 
this for the U (r) in leads to the relatively simple criterion for stability 


r? —6Mr + 8aV Mr — 3a > 0 (7.8) 


This is soluble analytically for r, but perhaps more illuminating to is to look at a 
plot of risco (the r corresponding to the innermost stable circular orbit) versus spin, 
as shown in Fig. Ifa = M, an ‘extremal’ black hole, then the solutions are 
r= M andr = 9M. Ifa = 0, we get the single solution r = 6M corresponding to 
the Schwarzschild case. The former of these shows that stable circular orbits (in the 
prograde direction) persist right up to the event horizon for a black hole spinning at 
the maximum mathematically allowed rate, i.e. a = M. We note, however, that various 
effects are likely to intervene before this maximum spin rate can be achieved (mainly 
the counteracting torque felt by the hole in absorbing radiation from the accretion 
disc [70]), and a = 0.998M is considered the most likely maximum attainable value, 
leading to a minimum attainable risco of 1.24M. 
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Radius of innermost stable orbit vs. spin 
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Figure 14. Plot of radius of innermost stable orbit vs. spin. 


More generally, by eliminating r between and the expression for k in 
we can get a useful relation between the black hole spin and the efficiency of energy 
release for a particle which has reached the ISCO. For prograde orbits this relation is 


a 2(2v2vI-# -k) 
M 3A- k) 


(7.9) 


This provides the plot of efficiency versus spin shown in Fig. Again taking the 
maximum attainable value of a/M as 0.998, this leads to a maximum attainable effi- 
ciency of 32%. 


7.2 Current Observations of Accreting Black Holes contd. 


The accreted matter can be supplied by a companion star, in the case of a BHB, or 
from surrounding gas in AGN. How much is captured by the black hole in the AGN 
case will depend on the density, temperature, magnetic field and angular momentum 
of the gas. The case when there is no angular momentum and magnetic fields can be 
ignored is known as Bondi accretion, which assumes that the gas can be treated as a 
fluid. Basically gas within the Bondi radius, defined as the radius beyond which the 
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Efficiency of energy release vs. spin 
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Figure 15. Plot of efficiency of energy release vs. spin (prograde orbits). 


gas particles have the escape velocity, will be captured and this can be many orders of 
magnitude larger then would occur if the particles were free particles with no interac- 
tions or collisions. Zero angular momentum is a rather extreme assumption and may 
not apply to any real system, but it presents an easily calculated order of magnitude 
estimate of the accretion rate. Angular momentum is likely to reduce any inflow rate 
and could even choke off the flow. Indeed accretion flows can be unstable when heat- 
ing by radiation from the central accreting mass is considered. Magnetic fields and 
turbulence also contribute to complexity and accretion flows are commonly observed 
to show chaotic variability. 


The accretion disc will be thin, with the ratio of thickness to radius, h/r < 0.1, 
when the luminosity is above about 1% Lgqa. At lower accretion rates the disc thick- 
ens as the collision rates drop which can then cause the efficiency to drop with further 
decreases in M. (The black hole swallows both the matter and energy released.) The 
properties of such radiatively inefficient flows (RIAFS) depend on whether the elec- 
trons (which radiate) and the protons (which have much of the energy) remain cou- 
pled, i.e. have the same temperature. There are still theoretical issues remaining as to 
whether this happens in such a highly magnetized plasma, but observations do indicate 
changes in behaviour in objects with low accretion rates. Not least is a tendency for 
jets to become common. The low luminosity of Sgr A* suggests it is in an extreme 
RIAF state. 
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The luminous accretion flows, however, are usually above 1% Legqa and so are thin 
and dense and optically thick resulting in much radiation being emitted as a black 
body. The spectrum of the disc is then produced by the sum of blackbodies from 
different radii which gives a v~!/3 spectrum at low energies up to the highest energies, 
which come from the smallest radii, where it turns over into a Wien tail. The peak 
temperature scales as Mrg '/4 This means that the most massive luminous black hole 
systems have spectral peaks in the UV, while stellar mass black holes peak in the X-ray 
band. 

Significant power can also emerge as a results of magnetic fields wound up in the 
accretion flow. They can emerge from the disc, become twisted and create a mess of 
reconnecting field lines powering a hot corona, with a temperature which can easily be 
100 times hotter then the disc itself. The million degree Solar corona with associated 
flares and outflows is a nearby example but even low luminosity AGN have a power 
10! times larger than the Sun, so the analogy should not be taken too literally. The 
corona is the source of most of the X-rays seen from AGN. Evidence discussed later 
indicates that the corona in bright objects is located above and close to the black hole 
(~ 10rg). The magnetic field structure also appears to power jets in some objects, and 
outflows in others. The degree to which the spin of the black hole is important (and 
maybe also a small ISCO) are still debated. 

BHB have stellar-mass black holes and accrete from their companion star. Disc for- 
mation is implicit due to the orbital angular momentum. A few sources such as Cyg 
X-1 are persistent, with their long-term luminosity varying by less than an order mag- 
nitude, but sometimes switching states from hard (coronal dominated) to soft (thermal 
disc dominated). Many of the BHB are however transient, having outbursts every few 
years or decades, during which the luminosity can rise from 10°? to 10°? erg s7! . Out- 
bursts follow a similar pattern in the luminosity-colour (hardness) plane, starting hard 
as the luminosity builds up towards Legga, then switching to being soft as the thermal 
disc emission dominates over the corona, before subsiding back to the hard state as 
they decline to quiesence [59]. At the peak of an outburst the nearest of these sources 
(eg A0620-00 at 2.5 kpc) can become the brightest in the sky, outshining Sco X-1 and 
even the Sun (except when it is flaring) by a large factor. 

There are also many similar X-ray binaries in our Galaxy which host a neutron 
star instead of a black hole. We know about the neutron star by observation of rapid 
pulsations due to the spin of the highly magnetic field, or from X-ray bursts which 
are thermonuclear flashes taking place very near the neutron star surface. Many of 
these are also transient and behave in a similar way to BHB. The depth of quiesence 
has been proposed as evidence of an event horizon in BHB by Ramesh Narayan and 
colleagues [47]. Black holes have no hard surface and can accrete matter and radiation 
whereas neutron stars must emit the gravitational energy released. BHB generally have 
lower quiesent luminosities than neutron star binaries which fits in with this picture. 
Whether these systems have similar enough accretion rates in quiescence for this to be 
a proof has yet to be determined. 
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Both BHB and some neutron star systems also show quasi-periodic oscillations, 
discovered by Michiel van der Klis in 1986 [74]. This phenomenon is complex and 
may be a rich source of information about the accretion flow and geometry, but it has 
yet to be fully understood and interpreted. To understand at least one of the frequencies 
we would expect to be seen, we now calculate the expected orbital frequency seen by 
a distant observer for material in the innermost orbit. 


7.3 Velocities and frequencies 


For astrophysical black holes, further important quantities for circular orbits are the 
velocity in the orbit, and the orbital period. For the former, the velocity determines 
(essentially via special relativistic effects) an important component of the shift in fre- 
quency which occurs for line radiation coming from the accretion disc, and for the 
latter, we can in some cases see directly the modulation of radiation associated with 
an orbital period. These are thought to provide one of the frequencies within the 
‘quasi-periodic oscillations’ which are occasionally observed from accreting stellar 
mass black holes and can have frequencies in the several hundred Hz region [74]. 

Working with a Schwarzschild black hole, then if we let v be the ‘ordinary’ velocity 
in a circular orbit, and a be the corresponding ‘rapidity’, then (temporarily dropping 
factors of c since this makes the formulae clearer), we have 


v=tanha and rọ = yv = sinha (7.10) 


We can thus use for h and the definition r?°ġ = h, to calculate v as a function 
ofr, obtaining the simple result 


2_ GM 

— r—-2GM 
For an object rotating in the innermost stable circular orbit (ISCO), at r = 6GM, we 
therefore find that its velocity is half the speed of light. (Note that, perhaps surpris- 
ingly, exactly the same velocity is obtained for the ISCO about a spinning black hole 
with any magnitude for its spin.) 

As regards orbital frequency, then since the t coordinate is the time as measured by 
a stationary observer at spatial infinity, if we can form do/dt then setting this equal 
to 27Vop will immediately give us the orbital frequency, in Hz, as measured by an 
external observer. 


ù (7.11) 


Thus using 
dọ _ dọdr_ $ ;_h GM 1 
E Ra =z 7.12 
dt drd È’ á r2 3 (1— 3GM) (7.12) 
and A l 
i= = (7.13) 
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do _, = ġọ_ [GM 
u Tob = 7 = a 


which despite being derived in full GR is exactly the Newtonian (Keplerian) result. 

As an example, evaluating at the innermost stable circular orbit, r = 6GM/c’, for 
a 10 solar mass black hole, yields the impressively high result of Vow = 218 Hz. Thus 
an emitting blob orbiting at this radius could potentially produce this frequency in the 
frequency spectrum of the source. 


we find 


7.4 Further AGN properties 


AGN range in luminosity from below 10!° to 10'+Lo and are situated at the centre 
of the host galaxy. They appear to respect the Eddington limit and show strong evo- 
lution in luminosity with redshift (and thus cosmic time). The quasar and therefore 
mass-building phase peaked at about redshift 2.5 with quasars being observed back to 
redshift 7.1. An important issue with AGN is that a plentiful supply of fuel leads to 
obscuration of the source itself. If the absorbing column density is Thomson thick (ie 
most photons are scattered by electrons on passing out of the source), then only heavily 
reprocessed emission is all that remains, emerging predominantly in the infrared. 

Type I AGN are largely unobscured and Type II obscured. The difference may be 
geometrical and due to our line of sight. In a broad class of source the obscuring 
medium is in the form of a thick torus at a radius of about a parsec. If viewed in the 
plane of the torus we see a Type II obscured object and if we view it face-on then it is 
Type I. If at the same time there is a jet, then face-on down the jet the object appears 
highly beamed as a blazar [72]. 


8 Measurements of the Masses of Black Holes 


The mass of a compact object can be measured from the motions of stars or gas in its 
vicinity using Newton’s Law of Gravitation. If a star makes a circular orbit around 
a massive black hole with velocity v at radius r then the mass of the black hole 
M = v’r/G. Multiplying factors emerge for different shaped orbits but the essence 
of measuring velocity and radius remains the same. For the black hole at the centre of 
our Galaxy, Sgr A*, the 3D motion of many nearby stars has been mapped to give the 
best estimate of 4 x 10°Mo. 

With stellar mass black holes orbiting a normal star, then constraints on the mass 
are obtained by measuring periodic doppler shifts of spectral lines of that star to yield 
its radial velocity curve of semi-amplitude K and period P. These are combined to 
give the mass function 

(Mp sini)? — PK? 
(M; + Mion)? 27G 


f= (8.1) 
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Further progress requires estimates of the mass of the normal star, perhaps from its 
spectral type, and of the inclination of the orbit 7. If the system is eclipsing then i 
must be high, but if not then other means must be used such as observation of periodic 
ellipsoidal variations in the lightcurve of the star due to tidal distortions. These can 
be modelled to yield the inclination. There are now several dozen black holes with 
masses estimated to about 10% in this manner. 

The mass of supermassive black holes (SMBH) is obtained from the motions of 
nearby stars and gas if they lie in relatively nearby galaxies where the region of in- 
fluence of the black hole can be resolved. This is roughly the radius within which 
the black hole-induced motions are faster than the general motions in the galaxy. The 
masses of more distant active SMBH can be determined using the reverberation tech- 
nique. Delayed variations in the ionizing UV radiation from the nucleus are seen in or- 
biting broad-line clouds — the whole system appears to reverberate with delays rang- 
ing from days to weeks. (The clouds show doppler-broadened emission lines due to 
orbiting close to the SMBH.) The light travel-time delay then gives the radius and the 
line width of the clouds gives the velocity required to determine the central mass. As- 
sumptions need to made about the geometry but careful work of the past two decades 
has established the success of the method |55]. It has enabled correlations to be built 
between the width of selected emission lines, optical luminosity and black hole mass. 

Current results indicate a maximum mass for observed neutron stars at just over 
2Mo, while stellar mass black holes have masses from ~ 4 — 20Mo, with an average 
of about 7Mo. Observed supermassive black holes have masses ranging from just 
below 10° to 2 x 10!°M,. It is possible that some UltraLuminous X-ray sources 
(ULXs), which have luminosities greater than the Eddington limit for a 10M. black 
hole, are Intermediate Mass Black Holes (IMBHs) with masses of 10? — 10° Mo. 


9 Measurements of Black Hole Spin 


Whilst mass measurements can be made at a large distance from an object, spin re- 
quires a probe which is close in, within 10r 4] The standard approach is to identify 
the inner edge of the accretion disc with the ISCO and then convert that radius to spin. 
This means that the inner disc must be luminous and detectable. 

Bright Galactic X-ray binaries in the soft state generally have a quasi-blackbody 
X-ray spectrum which can be used to measure the emitting area using the Stefan- 
Boltzmann law in an analogous way to that by which stellar radii are determined 
(blackbody luminosity L = 47. R?oT*) [B9]. Unlike the surface of a star, an accre- 
tion disc is not all at the same temperature, so the emission profile of the disc must 
be assumed and the black hole mass and distance determined in order to measure the 
ISCO and thus spin. 


4 Tf the mass accretion rate can be measured (from observations of the outer disc, say) and the total 
luminosity is known, then in principle spin can be deduced from the radiative efficiency using the 


relation in Fig [15] 
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Light bending effects are very significant for the accretion disc near a black hole [43], 
so we now give some details of light bending effects in GR, first in the Schwarzschild 
metric, and then for Kerr black holes. An impression of what we would actually see 
for the accretion disc around a Schwarzschild black hole, if sufficient resolution was 
available, is shown in Fig. [16] 


Figure 16. Pseudoimage of a luminous thin accretion disc around a black hole with 
doppler and gravitational redshift effects included. Light bending causes the distant side 
of the disc to appear as an arch above the hole. The smaller loop near the centre is due to 
scattering of photons from the (unstable) photon orbit at r = 3G'M/c’. The colour scale 
indicates flux received. 


This can be compared with the more distorted and asymmetric image we would see 
from a black hole rotating at the expected maximum rate of a = 0.998, in Fig. [17] 


Figure 17. As for Fig. [16] but for a Kerr black hole rotating at the maximum attainable 
value of a = 0.998M (see Section FI. The unstable photon orbit is now near r = 
GM/c’. 


9.1 Equations for photon motion and redshift 


The other results we need for understanding astrophysical aspect of black holes, are 
the equations of motion for photons, and results for how their energy changes during 
propagation. 
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For the equations of motion, we can use the above analysis, but with the interval 
function G(a", <") set equal to 0 rather than 1, once its functional variation has been 
taken. 

Tracking through the changes this causes, we find equation is replaced by 


2 
i = ke k (1 c) (9.1) 


r2 re 


and for the ‘shape’ equation, (3.9), one finds that it is now the ‘Newtonian’ term that 
disappears on the r.h.s., and we just get 

u 3GM > 

= ur. 


dp "TTZ 


(9.2) 


One can immediately confirm from this equation that there is a circular photon orbit 
atr = 3GM/c’, as mentioned above, but is it stable? 
To do the stability analysis, we rewrite the energy equation as 


t% 1 
T+ Venlr) = 5, (9.3) 
where b = h/ck, u = GM/c? and the effective potential 
1 2 
Vaar(r) = -3 (1 = 3 (9.4) 


Let us look at a plot of this function, Fig. We can see Vere(r) has a single maximum 
atr = 3u where the value of the potential is 1/(27*). This shows that the circular 
orbit at r = 3u is unstable and in fact no stable circular photon orbit is possible in 
the Schwarzschild geometry. 

Another immediate use of equation is in connection with light bending. Refer- 
ring to Fig. we can see that a suitable first solution in which the term 3G'M wu? /c? 
is completely ignored, is 

sind 
u=- p (9.5) 


where R is the radius of the body the gravitational deflection due to which we wish to 
work out. We iterate this equation by putting sin? p/P? for u? on the rhs. of (9.2), i.e. 


Cu 3GM 


igi! = am sin? ġ. (9.6) 


This is satisfied by the particular integral 


3GM 1 
ui = 22 R2 (1 + — COS 20) 5 (9.7) 


Astrophysical Black Holes 39 


Effective potential for photon motion 
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Figure 18. Plot of the effective potential for photon motion in the Schwarzschild geome- 
try. (Units of Veg are 1/(27p°) where p = GM/c?.) 


Figure 19. Gravitational deflection of light 


and adding this into the original solution yields 


sind  3GM (1 


1 
R | 22m? £0820) . (9.8) 


3 
Now consider the limit r — ov, i.e. u — 0. Clearly we can take sing = ọ, 
cos 2 ~ 1 there, and we obtain 6 = —2GM/(c’R) so that the total deflection (see 

figure) is 
4GM 

A¢g = =. 9.9 
o= R 0.9) 
This is the famous gravitational deflection formula. For the Sun it yields 1.77 seconds 
of arc, and was first verified in the 1919 eclipse expedition. More recent high precision 
tests use radio sources, since these can be observed near the Sun, even when there is 
no eclipse, and there is now no doubt that the GR prediction (which incidentally is 
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twice what had previously been worked out using a Newtonian approach) is accurate 
to a fraction of a percent. 

In the black hole context, much stronger bending can occur, for which this first 
approximation is insufficient, and the full equations of motion need to be integrated. 
This can lead to some interesting effects, illustrated below in the context of a Kerr 
black hole. 

The other important aspect to discuss is gravitational redshift. Quite generally, for 
a metric with a Killing vector field g+, the redshift due to gravitational time dilation 
will be given by the ratio of p: gą at emission and reception, where p is the photon 
4-momentum. g, becoming null, i.e. gu going to zero, marks the point where infinite 
redshift is possible since p can then be in the direction of q,, e.g. for counter-rotating 
Photons at the Kerr Stationary limit. This redshift from p: gı acts multiplicatively on 
any further redshifts which are occurring at a point due to special relativistic effects. 
E.g. for a point in an accretion disc, there will be a doppler redshift due to its motion 
relative to the initial path of an emitted photon, and then this will be multiplied by 
the gravitational redshift the photon experiences between that point and infinity. The 
combined effects of the redshifts can be seen in pattern of received flux seen in Figs. 
and\I7\ which are for the case where the accretion disc itself has constant emissivity. 


9.2 Light bending around a Kerr black hole 


The treatment of photon orbits for the Kerr case can proceed in the same way as in 
Section|7.1] which dealt with massive particles, except that now the interval function 
G(x”, &") evaluates to 0 instead of 1. Several textbooks, e.g. and [22], deal with 
the details of the resulting orbits, and in particular show that as in the Schwarzschild 
case, there are still no stable circular photon orbits. 

The possible orbits for non-equatorial photons can be quite exotic. As an example, 
we show in Fig. [20|the path of a photon emitted from an accretion disc in the equatorial 
plane, moving towards the black hole but with an initial upward component to its 
velocity. This moves round behind the black hole and is eventually able to escape to 
infinity, moving in a roughly opposite direction to its initial motion. This example is 
for a Kerr black hole with the maximum expected value for a of 0.998M — the ISCO 
and horizon radii are also shown. 


9.3 Iron Line Emission 


The method discussed above for Galactic X-ray binaries does not work for AGN since 
the much higher black hole mass means that the blackbody disc emission emerges in 
the far UV band which is unobservable because of absorption by Galactic hydrogen 
clouds. Observations are made instead of broad iron-line emission produced by X-ray 
reflection. This is the fluorescent and back-scattered emission produced by irradiation 
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Figure 20. Path of a photon emitted from an accretion disk that eventually escapes from 
the black hole in a direction about 180° from its initial direction. 
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Figure 21. Schematic cross-section of an accretion disc around a black hole with a central 
corona providing irradiation. The production of the reflection spectrum is illustrated. 


of the disc by the coronal X-ray continuum — see Fig.21] X-rays absorbed in the disc 
lead to line emission, particularly iron K emission at 6.4 to 6.95 keV, depending on 
the ionization state of the emitter. At higher energies the X-rays are more likely to be 
Compton-scattered back out of the disc rather than absorbed, producing what is called 
the Compton hump in the reflection spectrum. 

The spectral features are produced in the innermost parts of the accretion disc where 
gravitational effects are strong and lead to large doppler shifts and gravitational red- 
shifts which skew and broaden emission lines — an example is shown in 
Fig. Detailed modelling of the line shape, and generally the shape of the whole 
observed reflection spectrum, leads to an estimate of the inner radius of the disc and 
thus the spin. Radii are determined in this method directly in gravitational units so the 
distance and mass of the black hole are not required. 

The broad iron line method can also be used in certain states of BHB and has 
been found to give consistent results with the blackbody method. 

Confirmation that the reflection modelling is correct has recently been obtained from 
X-ray reverberation [73]. This involves the small (tens of seconds to tens of minutes) 
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Figure 22. Upper: The unblurred reflection spectrum including the iron-K emission line 
and Compton hump (solid line) appears relativistically blurred to the observer as shown 
in the line fitted to the data points. The obvious dominant redshift is gravitational due to 
the proximity to the black hole. The source is the Seyfert galaxy Mkn335, observed by 
NuSTAR during a low state. Modelling of the data shows that the black hole in Mkn335 
has a spin a/M > 0.97. Most of the emission seen in the spectrum emerges from within 
2rg of the event horizon [53]. Lower: Time lags measured by NuSTAR from the AGN 
SWIFTJ2127 showing both iron K and Compton peaks. 
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time delays or lags whch occur between the detection of an intrinsic change in luminos- 
ity of the corona and the corresponding change in the reflection. The energy spectrum 
of the high frequency time lags shows a broad iron line indicative of reflection (lower 
panel of Fig. [22). 

The results for a number of objects imply that high spin is common for the brightest 
observed sources and the corona in those objects is compact and at a height of 5— 10rg. 
In this regime, strong light bending occurs such that much of the coronal emission is 
bent down to the disc making the reflection strong. Month-long temporary episodes 
have been seen in a few sources in which the continuum drops to being almost unob- 
servable, with only slight changes to the reflected spectrum. It appears that the corona 
has dropped below 5r, in these cases and that light bending is preventing much of the 
continuum being observed at all. Most of the reflection spectrum is then emerging 
from within 2 — 3rg, which is the closest emission easily seen from around any black 
hole [53]. 

Further evidence of a compact corona is obtained from microlensing studies of dis- 
tant quasars which are lensed by intervening galaxies. This produces several images 
of the quasar which follow different light paths through the galaxy. Stars in the galaxy 
crossing one of those paths amplify the image by lensing to a degree depending on 
how the angular sizes of the emission region of the quasar and the Einstein ring of the 
star match one another. For typical values this means that higher amplitude variability 
is expected from smaller emission regions. Both optical and X-ray observations have 
been made of several lensed quasars which show that the X-ray emission region is 
much smaller than the optical one, and in several cases is less than 10r, in size [B]. 


10 Future Observations of Astrophysical Black Holes 


Many exciting developments are expected over the next decade or so in observational 
studies of black holes. One important step will be the networking together of several 
millimetre-wavelength telescopes to form the Event Horizon Telescope which will 
resolve the “shadow” of the black hole in Sgr A* [13]. Extreme light bending of 
the emission from the accretion flow onto the black hole should reveal a dark patch, 
the shadow, at the position of the black hole. The diameter of the shadow is determined 
by the paths of light rays emitted from just above the event horizon to the observer and 
corresponds to about 10r, almost independent of the black hole spin. 


Calculation of the shadow radius: 
We use the results above for the photon effective potential in the Schwarzschild case. 


Putting together equations and (9.4), we find 


2 
= : (1 x) — (10.1) 


h r r b2 
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where u = GM /c? and b has been defined as h/kc. Note the dot here is not derivative 
w.r.t. proper time, which is zero for a photon, but w.rt. an affine parameter along the 
path. We do not need to identify this parameter explicitly, since we are going to work 
with the shape of the orbit. Specifically, noting r?°ġ = h, we have 


2 1/2 1/2 
ii — 1 dr = 1 1 1 2u (10.2) 
r4 r? dọ B2 r? r 


Solving (10.2) for u = 0, which is the limit where the only role of the central mass 
is to define the origin of the (r, ġ) system), we get r = b/ sin ọ. This is the equation 
of a straight line, which if it continued all the way past the central object, would pass 
it with perpendicular distance b, so we can identify this b as the impact parameter. 
This generalises what we did above, where we considered the case of a light ray just 
grazing the surface of an object of radius R, for which the unperturbed solution was 
given by (9.5), ier = R/sing. 

Now we consider the exact equation again, and relate to a black hole. We know 
from the plot of photon effective potential in Fig. that photons straying inside r = 
3GM/ c are not free to escape to infinity. We therefore want to use this distance as the 
radius of an ‘object’ which the incoming photon just grazes. Using R for this object 
radius, we have 


This tells us that 


1 1 nN R N? 
=—(1 e b= 
4 val 2) K R( az) 


This expression is valid generally, and e.g. tells us that the Sun’s effective shadow is 
about 3 km bigger than its coordinate radius, but applied with R = 3GM/c for a 
black hole, we find b = 3V/3GM /@, which is therefore the radius of the geometric 
shadow cast. Interestingly, the main effect of including black hole spin in the calcula- 
tion is to cause a shift in the centroid of the shadow rather than its size — see e.g. 
for details. 


For Sgr A* this distance we have just worked out corresponds to an angular diameter 
of 52 microarcsec. This is 2.5 x 10~!° radians, and is very small indeed. Nevertheless, 
the angular resolution, approximately equal to the ratio of observing wavelength to 
telescope diameter, is about 20 microarcsec for intercontinental mm VLBI where the 
effective telescope diameter is 10,000 km. 

Interestingly, the 6 billion Mo black hole in M87 is 1500 times more massive than 
Sgr A* and 2000 times further away which makes its event horizon appear only slightly 
smaller, at 40 microarcsec. Plans are also underway to resolve the M87 black hole 
shadow, together with the base of its jet which should help us understand how it is 
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accelerated. All other black holes subtend smaller angular sizes. 3C273, for example, 
is 34 times more distant than M87, so its shadow will be about | microarcsec in size. 

Near infrared interferometry at ESO’s Very Large Telescope in Chile will use wave- 
lengths around 2 microns and effective diameters of a 100 m to give a resolution of 4 
milliarcsec. This will not resolve the black hole shadow but will resolve the motion of 
the orbiting stars much more precisely, particularly as the positional accuracy can be 
more than an order of magnitude better than the resolution. General relativistic effects 
on the orbits will be seen and fainter, possibly closer, stars resolved. 

Of great interest would be a pulsar orbiting Sgr A*, which if a millisecond pulsar 
could prove to be a very stable and reliable orbiting clock. None have been found yet 
and the ionized plasma that bathes the inner region is a problem for sensitive radio 
searches. A type of pulsar called a magnetar was found in 2013 within 3 arsec of Sgr 
A*, but this magnetar unfortunately appears to be rather unreliable as a clock [10]. 

A gas cloud, G2, was discovered on a highly eccentric orbit close to Sgr A* in 
2012 [56]. At pericentre in 2014, about 3 Earth masses of material passes within about 
2000r s. Ram pressure from the accretion flow should cause it to spiral inward over the 
following few years leading to an increased accretion rate and associated display. It is 
in a similar path to an earlier discovered cloud G1 and other clouds. Their apocentre 
and orbital plane lie close to the orbits of some of the surrounding luminous stars. The 
winds of the stars may be the origin of the inspiralling gas clouds. 

Further past activity of Sgr A* has been revealed by X-ray fluorescent emission 
lines seen from molecular clouds many tens to hundreds of light years out from the 
Galactic Centre. One interpretation suggests that the luminosity of Sgr A* may have 
been almost a million times higher about 100 yr ago [7]. 

Sufficient numbers of black hole spin measurements are expected from X-ray stud- 
ies of AGN that the spin history can begin to be understood. Present work indicates 
many high spin objects in the mass range from 10° — 5 x 10’ Mo. The accretion flow 
within 10r, and particularly within 3r, will be mapped and the corona understood. 
Indeed the whole power output of accreting black holes in terms of radiation, winds 
and jets will be explored. 

Tidal disruption events due to stars straying close to a massive black hole are begin- 
ning to be observed. Black holes with masses above 10° Mo can swallow stars whole 
in the sense that they are not tidally disrupted until within the event horizon. The 
tidal forces of the more common lower-mass black holes can however destroy stars at 
greater distances leading to the production of a short lived, very luminous, accretion 
disc which decays with a timescale of about a year. 

Studies of AGN feedback will accelerate with use of telescopes such as the At- 
acama Large Millimtre Array ALMA, the James Webb Space Telescope JWST, the 
Large Synoptic Survey Telescope LSST and the Advanced Telescope for High Energy 
Astronomy Athena. The growth of massive black holes through accretion and mergers 
will become understood. 

Actually testing GR with black holes is difficult. The astrophysical phenomena 
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observed so far are explained within the Kerr metric. Some precision tests may occur 
if a pulsar is serendipitously found in a close orbit around a stellar mass black hole 
in our Galaxy or around Sgr A*. The discovery of gravitational waves from merging 
black holes will of course transform this aspect of the study of black holes. 

More generally, the detection of gravitational waves emitted in events involving 
black holes has great potential for tying down aspects of astrophysical black holes that 
are poorly understood at present. These include the numbers of intermediate mass 
holes in the range between stellar mass black holes and SMBHs in galactic centres, 
and the evolution of the number density of black holes with cosmic epoch. Both these 
aspects are key targets for a future space-bome gravitational wave detector mission, 
such as eLISA or a future variant, and form a key component of its science mission 
[17]. There are many further aspects to the importance of gravitational waves in black 
hole astrophysics, but we defer coverage of these to other chapters of this book. 


11 More general spherically symmetric black holes 


Real black holes are embedded in a universe in which we know there are effects match- 
ing those of a cosmological constant, accounting for about 0.7 of the total energy den- 
sity. Thus it is of interest to embed a black hole in a de Sitter universe. Additionally, 
black holes in anti-de Sitter universes are very interesting objects theoretically, since 
via the AdS/CFT correspondence (see e.g. ), the strong gravity effects of black 
holes can be replaced by weak perturbative effects in a conformal field theory living 
on the boundary of an anti-de Sitter space (though we note that the dimensionality 
involved does not necessarily match that of 4d spacetime). 
Returning to the derivations of Section we find that to satisfy the new field 
equation 
G uv + AGuy = —8aT yy (11.1) 


arising from introducing a cosmological constant A, the new equations the Ricci tensor 
has to satisfy are 
Ry = Agu (11.2) 


Using the components already given for the Ricci in and above, one finds that 
this leads to a simple modification to the Schwarzschild values of A and B. We find 
that B is still equal to the inverse of A, and that now 


-|r (11.3) 


These are therefore the metric coefficients for a black hole of mass M embedded in a 
de Sitter universe, if A is positive, and anti-de Sitter, if A is negative. 

Horizons occur (for our current form of the metric) when B — co, hence here 
where A = 0. IfA > 0, this now leads to an additional horizon being present, beyond 
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the Schwarzschild one near the black hole centre, at a position close to r = \/3/A (we 
are adopting a convention here in which A has dimensions length~*). This is a version 
of the ‘de Sitter horizon’ — its exact position is modified slightly by the presence of the 
mass. 

A further generalisation we can make, is to the case where the black hole is charged. 
As mentioned earlier, this is less compelling physically, since we do not expect signif- 
icant charge separation to have occurred in the formation of a black hole, but the 
solution for this type of black hole, called the Reissner-Nordstrom solution (see e.g. 
for references), is still of great interest theoretically. This also falls into our gen- 
eral scheme, and in fact the most general spherically symmetric black hole, which we 
can dub RNdS, for Reissner-Nordstrom-de Sitter, is given by 


2GM Ay ¢ 

re 3 Tp 
where q is the black hole charge. Note there are now (in general) three horizons. The 
Reissner-Nordstrom solution on its own introduces two horizons — one like the normal 
Schwarzschild one, and another usually much smaller one associated with the charge 
— and then for A > 0 there is an additional de Sitter horizon. 

An advantage of dealing with all these cases via the A, B form of the metric, is that 
we can give a unified treatment of energies and angular momenta in orbits. Repeating 
the analysis of Section Bl for this new case, we already have ù from and can get ï 
from 


B=A'!, A=1 


(11.4) 


ee ee a ee 
dr 0 dt i i a 
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For a circular orbit, both ¢ and ¥ have to vanish and this combination of expressions 
therefore says that 


where A' = dA/dr. Combining again with the expression for * = 0 we can solve for 
k? and h?, obtaining 
2A? PA! 
k = ——_,_ and h? = ———_ 
2A—rA’ 2A—rA’ 
So we now have the energy and angular momentum for a particle in a circular orbit 
for quite a wide range of metrics. To calculate criteria for stability, one can show that 
it is useful to work in terms of the quantity (see Lasenby (2014) in preparation) 


2 
T(r) = (=) = tanh? u (11.5) 
c 
where v is the velocity in a circular orbit at radius r and u is the corresponding rapidity 
parameter. From what we have already derived, one can show 
rA' 


T= 57 and rọ = yv = sinhu (11.6) 
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and then the criterion for stability turns out to be that 
2T? — rT' = 2T (11.7) 


should be negative. This simple criterion applies across all the spherically symmetric 
black hole cases. For example, inserting the A appropriate for a Schwarzschild-de 
Sitter metric, , we find the stability criterion 


18M? — 15MAr? —3Mr +4Ar* <0 (11.8) 


for stable circular orbits (in units with c = G = 1). For appropriate astrophysical 
values of M and A, and restoring units, the outer solution to this is well approximated 


by 
3GM \ 13 
Tstab © G ) (11.9) 


We can apply this formula for any central mass, not just a black hole, and one can 
show that e.g. for a cluster of galaxies with M = 105 Mọ, and for the measured 
value of the cosmological constant, then Tsab ~ 7.2 Mpc. This is interestingly close 
to the maximum cluster sizes observed, though of course circular orbit stability is not 
likely to be the direct criterion one would apply in this case, where effects on velocity 
dispersion would be more appropriate (see for further details). 


12 Primordial black holes 


Primordial black holes could have formed as a result of the high densities present soon 
after the Big Bang [8]. Their masses would be comparable to the particle horizon mass 
at their formation so it can range from the Planck mass (1078 kg) at the Planck time 
(10-*% s), to 10° Mz at a time of 1s. Primordial black holes of mass 10!? kg would 
have formed at 107% s and would be evaporating now due to the emission of Hawk- 
ing radiation (see below), which would peak at 100s of MeV. Observations provide 
a limit on the background intensity of 100 MeV gamma-rays such that evaporating 
primordial black holes cannot account for more than 1078 of the critical density of 
the Universe [52]. The frequency, energies and timescales of the final explosions do 
not match typical observed gamma-ray bursts. It is therefore unlikely that they are an 
important constituent of the Universe. 


12.1 Hawking radiation 


An initially very surprising result about black holes was proved by Hawking in 1974 
. This is that despite their name, the horizon radiates energy, as though it was a 
blackbody at the temperature of 


he 


Tean — = 12.1 
8rkgeGM paN 
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where kp is the Boltzmann constant and M is the black hole’s mass. The types of 
black hole we have mainly considered have a minimum mass of several Mo, and so 
their temperatures are less than ~ 6 x 107? K, and therefore negligibly tiny. However, 
a 10! kg primordial black hole of the type mentioned above, would have a temperature 
of ~ 1.2 x 10!! K, and would therefore produce highly energetic radiation. 

The emission process as calculated originally by Hawking was explained in terms 
of particle production in a second quantised treatment of a massless scalar field sur- 
rounding the black hole, but the principles extend to all kinds of field, massive or 
massless, which will be radiated with ‘statistics’ (i.e. Bose-Einstein or Fermi-Dirac) 
as appropriate to their spin, and with a thermal distribution with the temperature 

12.1). 

A heuristic way of demonstrating this effect is to consider vacuum fluctuations in 
the space just outside the horizon, which will lead to the production of virtual parti- 
cle/antiparticle pairs there. Sufficiently close to the horizon, one of the pair will be able 
to ‘tunnel’ through the horizon and inside can become a real particle with negative en- 
ergy. (The possibility of negative energy arises since as we saw above (Section Z1), 
the timelike vector in the t-direction tips over to become spacelike inside the horizon, 
which coincides with the Stationary Limit surface in the Schwarzschild case, and so 
the projection of the particle momentum onto the local t-vector may be negative.) The 
particle remaining outside, which has positive energy, can therefore be emitted off to 
infinity, whilst the absorption of a negative energy particle by the black hole, decreases 
the mass M, thus providing the overall energy for the emission. 

This type of derivation can be used to calculate the energy of the emitted particle, 
which turns out to be independent of exactly where (though assumed to be close to the 
horizon) the virtual pair is produced, and yields an energy E = kgT which is only a 
factor of 27 higher than the correct result from (72.7) — see Chapter 11 of {30 for 
details. 

A further possibility for a calculation avoiding the full rigours of quantum field 
theory (QFT), is to work with a first-quantised field surrounding the black hole. For 
example, a detailed derivation for a massive Dirac field is given in Section 8 of [B5]. 
This yields precisely the result together with the correct Fermi-Dirac Statistics 
for a spin-1/2 particle, and the calculation can be extended to electromagnetic and 
scalar particles, with the same results that would be obtained from QFT in each case 
(Lasenby et al unpublished). However, it is only in a QFT context that one can be 
properly certain of the physics involved, and of how to go about making unique choices 
for the branches of analytic functions that are involved, and so a QFT approach is still 
necessary to be fully confident of the results. 


12.2 Link with surface gravity 


What a full QFT approach leads to is a very general result that singles out the ‘surface 
gravity’ of a black hole (of any type) as the important quantity. The surface gravity 
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would appear to be infinite at the horizon, since no particle can escape the black hole’s 
pull there. However, the definition of surface gravity intended here is that appropriate 
to the force felt at the horizon, but which would be evaluated by an observer at infinity. 
In e.g. , Exercise 6.4, this force is shown to differ from that at the horizon by the 
redshift factor from the horizon to infinity, which is also infinite. Specifically, one finds 
for a Schwarzschild black hole that the local force necessary to keep an observer of 
mass m stationary at radial distance r is (see e.g. Section 7 of [35]) 


GMm 2GM\~'/? 
F= 2 (1 va ) (12.2) 
The gravitational redshift factor between the horizon and infinity is (1 = ee 


hence multiplying by this, we obtain the simple expression GMm/r? as the force that 
the observer at infinity would ascribe as being necessary to keep the observer at radius 
r at rest. (In the exercise in Wald ( ), this difference is explained in terms of the 
different tensions necessary at the ends of a (massless) rope that extends from the 
observer at infinity to the one at radius r.) 

This may seem a mundane and wholly Newtonian result, and exactly what we would 
expect ‘at ininity’. However, we are now entitled to put r = 2GM/c’, i.e. the value of 
r at the horizon, into this expression, in order to evaluate the ‘surface gravity’ of the 
black hole. This is defined as acceleration per unit mass, and customarily denoted «K, 
so we have found 

A 


£= JGM 


The contribution of QFT at this point, is to tell us that a horizon with surface gravity 
k radiates particles with a thermal temperature 


(12.3) 


o Ak 
2rkge 


(12.4) 


Unruh in 1976 showed that this is the temperature of a thermal heat bath that 
develops around an observer accelerating at rate « in Minkowski space. We can now 
understand the link with surface gravity in black holes, at least heuristically. By the 
equivalence principle and working close to the horizon so that global curvature effects 
are not important, we can expect that the heat bath seen by a stationary observer feel- 
ing the effects of the black hole surface gravity «, should be the same as the heat bath 
seen by an observer in Minkowski space accelerating at the same rate, and we can see 
indeed that equations and agree with for this Schwarzschild case. 
The advantage of this route, however, is that if we work out the surface gravity in more 
complicated cases, such as for Kerr or Kerr-Newman or even de Sitter spacetimes 
(which as we have seen, also possess a horizon), then we can make the same transi- 
tion to temperature as via . For example, for our general spherically symmetric 
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metric (2.7), it is easy to show that the surface gravity is 

A’ 
k= ——— 
2 V AB B-!=0 
(The point where 1/B = 0 is picked out as being the position of the event horizon.) 
Evaluating this for the Reissner-Nordstrom metric discussed above, and for clarity 
temporarily putting G = c = 1 etc., so 1/B = A=1—2M/r + ¢@’/r’, we find that 
the surface gravity at the outer horizon is 


(12.5) 


Elia (12.6) 
(w+ VP-P) 


and the Hawking temperature of the radiation emitted in this case is given by inserting 


this k into (72.4). 


KRN = 


12.3 Astrophysical aspects of black hole evaporation 


We can compute the rate at which a Schwarzschild black hole loses mass via the 
formula for the luminosity of a blackbody at temperature T, i.e. cT* where o = 
mk /(60h3c’) is the Stefan-Boltzmann constant, together with the area of the horizon, 
which is 4x(2G.M/c?)?. This yields 


M oh 
15360 rG? M? 


Solving this differential equation, we find that the cube of the mass declines linearly 
with time: 


Me = —power radiated = (12.7) 


Ch i 
5120 7G? 
where Mo is the initial mass and t is measured from formation. This gives an evapo- 
ration time in terms of Mo of 


M3(t) = Mẹ — (12.8) 


tevap = 2.66 x 10° MG yrs (12.9) 


If the current age of the universe is 13.8 billion years, this means only black holes with 
a mass less than 1.7 x 10!! kg will have had a chance to decay by now. Due to the 
M~? dependence in the expression for radiated power, (72.7), nearly all the energy 
emitted is confined to the very last moments of the black hole’s life. 


12.4 Black hole entropy 


The fact that black holes radiate like black bodies, suggest that they should have an 
entropy. If we identify Mc? as the hole’s energy, then the thermodynamic relation 


dU = TdS (12.10) 
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along with our identification of the BH temperature in (22.1), yields (first working in 
units with G = h = c = kg = 1 for clarity) 


dM x 8M =dS, sothat S=4ntM? = a (12.11) 


where A = 161M? is the black hole area. Putting back the units, we find that the 
black hole entropy in units of kp, is 


Da (12.12) 


where lp = \/hG/c3 is the Planck length and A = 4r (2GM/ 2} is the BH area. For 
an astrophysical black hole, this entropy is large, and as we have seen scales like M?. 
The result of this is that it is thought that supermassive black holes in the centres of 
galaxies dominate the entropy budget of the observable universe. Estimates of this are 
given in [IS], with the result that SMBH may contribute up to 7 orders of magnitude 
more entropy than stellar mass BHs, and approximately 15 orders of magnitude more 
than the next largest component, due to photons. 

This identification of entropy with 1/4 of the event horizon area expressed in units of 
D; is the starting point for a great deal of work connected with whether a microphysics 
approach to black hole entropy is possible. In this approach we would seek to obtain 
the same result by counting states, and then assigning S = kgln N, where N is the 
number of available microstates. A related question is what happens to the informa- 
tion about the matter and radiation from which the hole was originally formed. This 
information is screened beyond the horizon, and our lack of knowledge of the internal 
arrangements and compositions of this original material provides a satisfactory un- 
derstanding of why the black hole should have an entropy at all, as already proposed 
by Beckenstein [2], before Hawking’s discovery of black hole radiation. However, at 
the end stage of evaporation of the black hole, we no longer have a screen, and unless 
the information is somehow encoded in the phases of the emitted particles and fields 
composing the Hawking radiation, so that this is not random after all, then the in- 
formation has been irretrievably lost. This would seem to violate the basic principles 
of quantum mechanics, which demand unitarity of transformations between beginning 
and end states, and would not allow the disappearance of information in this way, 
given that there is no longer a casual horizon beyond which the information is hidden. 

These are deep and fundamental questions, which we do not venture further with, 
but are discussed elsewhere in this volume. We note that our discussion so far has been 
just in terms of a Schwarzschild black hole, but equivalences of the same nature work 
for all types of black hole, and are encoded in the ‘Laws of black hole thermodynam- 
ics’. We now briefly discuss these concentrating on their astrophysical implications, 
and in particular the possibility of extracting energy from rotating black holes. 
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12.5 Laws of black hole thermodynamics and the Penrose process 


The generalisation of (12.10) to an uncharged rotating black hole reads (in natural 
units again, for clarity): 


dU = TdS + QydJ 


k (12.13) 
ie. dM = 5 -dA+ QudJ 
TT 


where Qq is the black hole angular velocity at the horizon, and we have again made 
the identification of temperature with surface gravity/27 and of entropy with one quar- 
ter of the horizon area. The new feature here is that work can be done to change the 
angular momentum J of the black hole. (Recall that this is related to the mass and 
spin parameters via J = aM.) 

Equation is known as the ‘First law of black hole thermodynamics’. There 
are equivalent versions for black holes of the other three laws of thermodynamics 
as well — see e.g. Section 9.9 of for a discussion of these. For astrophysical 
processes, the other law of immediate interest is the equivalent of the Second law, 
which states that in any classical process, the area of a black hole horizon, which we 
know measures the black hole’s entropy, does not decrease. The restriction to classical 
processes is necessary, since as we have seen, Hawking radiation succeeds in reducing 
both mass and surface area. 

The most interesting immediate application of these laws is to (classical) processes 
in which we attempt to remove energy from the hole. This cannot be done for a non- 
rotating hole, but the QydJ term in the First law provides a route through to this for 
rotating holes. 

The first process of this kind to be discussed, was the Penrose process [54]. In 
this process, an incoming particle enters the ergosphere of the black hole. It then 
decays into two particles. Particles can escape from the ergosphere (which lies outside 
the horizon), and we arrange the initial trajectory and decay such that one of the 
decay products escapes, and the other falls into the black hole. The key observation 
is now that inside the ergosphere the g, Killing vector corresponding to invariance 
of the metric under time displacements, becomes spacelike, (remember the definition 
of the ergosphere is where g? changes sign), and it is possible for a particle to have 
negative energy when its 4-momentum is projected onto it. If this happens, then from 
conservation of energy, the energy of the emitted particle will be greater than that of 
the original particle, and we will have succeeded in extracting energy from the black 
hole. This energy has come via the BH’s absorption of the infalling negative energy 
particle. 

One can also analyse what happens to angular momentum in this process, by con- 
sidering the projections onto the other Killing vector (corresponding to ġ invariance 
of the metric) gg. This reveals (see or [22\ for details) that the hole absorbs 
negative angular momentum as well as negative energy (meaning the emitted particle 
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enjoys a boost to its angular momentum as well as energy, relative to the incoming 
one) and that the black hole parameter changes obey the inequality 


6M > Qy6J (12.14) 


Now, the expressions in the Kerr case for the quantities appearing in the First law 


aie 
Horizon area. A=8n (m? + /Mt— P) 
J 
2M (M? + VM P?) (12.15) 
VIZE 
2M (m? $ VW T) 


(again see or for details). We thus have that the response of the horizon area 
to a general change in M and J is 


Horizon angular velocity Qy = 


Surface gravity K= 


OA OA 8rd 
= amo " aJ oJ = Qu VM? — J? (M —Qy6J) (12.16) 


From , we see that despite the fact that both M and J decrease, the r.h.s. here is 
always positive in the Penrose process. This is an example of the Second law of black 
hole thermodynamics in action — for a classical process, the BH entropy, as measured 
by the horizon area, must always increase. 


12.6 Adiabatic (reversible) changes 


An interesting case to consider astrophysically, is the limit of gradual changes in mass 
and angular momentum, which satisfy the thermodynamic notion of adiabaticity, in 
particular the changes are reversible. This is possible, for example, for particular 
versions of the Penrose process. In the black hole context, since horizon area equates 
to entropy, an adiabatic process will have 5A = 0. 
Solving (12.16) for 6A = 0, immediately yields the following differential equation 

for M(J): 

dM _ J 

dy" om (2 + V= F) 


for which the following (implicit) solution works: 


2M6 = M?(J)+4/M4(J) — J? (12.18) 


where Mo is a constant. As an example, suppose we start with a Schwarzschild black 
hole with mass Mo, and gradually feed in material with positive angular momentum. 


(12.17) 
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We will be able to do this until an extremal black hole with J = M? is reached, at 
which point the mass is V2Mọ. An interesting point is what happens to the horizon 
radius Touter in this process. We know the general expression for it (in terms of M and 


J =aM) is 
1 

outer = M 4 VM (m? + / M4 J?) (12.19) 
and naively we might expect this to increase as the hole gains mass. But in fact we can 
see for the adiabatic process just discussed, 


© A  16rM 2M? 
"outer = SM 87M M 


(12.20) 


The horizon radius therefore goes down, reaching a minimum of V2 Mọ at the end 
of the process (corresponding to the rower = M value for an extremal black hole), 
having started at 2 Mọ. 


12.7 Other processes for extracting energy from a spinning black hole 


As well as the Penrose process, there are other ways of extracting energy from a ro- 
tating black hole. One of these is ‘superradiance’, which is the analogue for waves of 
what happens for particles in the Penrose process. The possibility of such a process 
was first drawn attention to specifically in the context of the Kerr solution by Starobin- 
ski in 1973 Ñ68]]. It involves an incoming radiation field of the form 


p ~ polr, Oje em? (12.21) 


for which part of the wave (the ‘transmitted’ wave) is absorbed by the black hole, 
and the other part (the ‘reflected’ wave) reaches infinity again. In the same way as 
for the Penrose process, due to the spacelike nature of the g, Killing vector inside the 
ergosphere, the transmitted wave can have negative energy, so that the reflected wave 
carries increased energy to infinity (specifically it will have greater amplitude at the 
same frequency as compared to the incident wave). We can derive the condition for 
this to occur (at least in a heuristic way), as follows (see e.g. Section 8.8 of [22]). 

We can think of the wave as composed of quanta with energy hw and angular mo- 
mentum hm. When these are absorbed by the black hole, the change in the black hole’s 
parameters will satisfy 

oJ hm 

ôM hw 
But we also know that 6M — Qy6J > 0 from the fact the horizon area must not 
decrease. Putting these together leads to 


(12.22) 


ôM 
— (w—-Qym)> 0 (12.23) 
wW 
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and so the condition for superradiance (65M < 0) is 
0 <w < mAy (12.24) 


Further details of this process, which only works for bosonic fields, can be found in 
Section 12.4 of [75]. 

Our final example, concerns what may be an important way in which spinning black 
holes return energy to the environment, and which was referred to briefly above in 
Section|6.1] namely the Blandford-Znajek effect [4]. 

Here, one considers a Kerr black hole immersed in an ambient magnetic field, typ- 
ically associated with an accretion disc around the black hole. A possible picture of 
the process (see e.g. Section 8.9) represents the black hole horizon as a moving 
conductor within the magnetic field. Indeed the Kerr black hole horizon, in common 
with all stationary event horizons, can be modelled as having a resistance of 4m in ge- 
ometrical units (377 Ohms in ordinary units). This rotating conductor coupled to the 
magnetic field, generates a current which flows between the poles and equator — it is 
in this sense like a dynamo. The main analytic field configuration chosen by Blandford 
& Znajek was a ‘split monopole’ — i.e. a different sign magnetic monopole solution 
in each hemisphere, with magnetic field lines pointing radially outwards or inwards. 
This needs a current sheet on the join between them in the equitorial plane, which is 
presumed to be supplied at least outside the hole by currents in the accretion disc. 

The radiated power can be calculated from the Poynting vector, and gives similar 
results as for a rotating magnetic dipole pulsar model evaluated on the light cylinder. 
Specifically we find (see e.g. Section 14.3) 


24 02 2.4 02 
B BrrpQy B Bry Q 


Cc C 


By \ pax / M V 
~ 4.1 x 10” ergs! (=) 
ee (55) M) \10°Mo 


Here Bpr and By are the magnetic field values at the ergosphere and horizon respec- 
tively, rg and ry the corresponding r values, and Q ș the magnetic field rotation rate, 
which is typically ~ 5Qu. This EM power will presumably be manifested in a mag- 
netised particle wind and jet, removing angular momentum and energy from the black 
hole in the process. Subsequent work (see e.g. Komissarov ) has shown that the 
Blandford-Znajek (BZ) mechanism is stable, and clarified that important components 
of it take place within the ergosphere itself, emphasising its links with the Penrose 
process and superradiance. The BZ mechanism is now thought to be an important 
component in generating at least some of the high energy jets we see in astrophysics 


pEM 


(12.25) 
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13 Conclusions 


Black holes are now an integral features of our cosmic landscape. Many millions of 
stellar mass black holes reside in our Galaxy. A supermassive black hole lies at the cen- 
tre of all massive galaxies. We know little about most of them, unless they are massive 
and nearby or accreting gas from a stellar companion or the surrounding interstellar 
medium. The inner part of the accretion flow can be extremely luminous. making the 
immediate surroundings of the blackest parts of the Universe into the brightest. 

The behaviour of observed black holes can so far be explained well by General 
Relativity and the Kerr metric. The physics of accretion onto black holes, and of the 
outflows and jets that often accompany inflows, is complex and leads to complicated 
observational phenomena which lie at the forefront of astrophysics. Understanding 
how quasars work remains a significant astrophysical challenge. The consequences for 
a galaxy hosting a supermassive black hole are profound. It is likely that the energy 
released by the growth of the black hole plays a decisive factor in its final stellar mass 
and possibly its physical size. 

Black holes are intrinsically relativistic objects which have stimulated and tested 
physical understanding to the extreme. The internal structure of the black hole, within 
the event horizon, is beyond direct observation and we have not discussed it here, 
although it is clearly of great theoretical interest. In this Chapter we have looked at 
black holes as observable physical objects and considered how they work as engines 
of gravitational energy release. 

The future of research into astrophysical black holes is very bright, as more tele- 
scopes probe ever deeper into more wavebands to uncover new objects, features and 
phenomena. Serendipitous discovery plays a key role in the history of astronomy and 
we look forward to the discoveries of a millisecond pulsar orbiting close to a massive 
black hole, a rogues gallery of black hole shadows, a range of quasars above redshift 
10, a complete X-ray spectral-timing deconstruction of the innermost accretion flow 
around a rapidly spinning black hole and, most of all, some new phenomena of black 
holes that we have not yet anticipated. 
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